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Abstract 

Applying the general framework for local zeta regularization proposed in 
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1 Introduction. 


In Part I of this series mmm we have considered in general local (and global) 
zeta regularization for a quantized neutral scalar field in a d-dimensional spatial 
domain hi, in the environment of (d + l)-dimensional Minkowski spacetime. 

In the present Part IV we consider a massless field confined within a d-dimensional 
rectangular domain = (0,ai) x ... x (0,0^), with Dirichlet boundary conditions; 
applying the general scheme of Part I we renormalize several observables, namely: 
the vacuum expectation value (VEV) of the stress-energy tensor, the pressure on 
boundary points, the total energy VEV and of the total force acting on any side of the 
box. All these renormalized observables are represented as sums of series converging 
with exponential speed, for which we give quantitative remainder estimates. Our 
results hold for an arbitrary spatial dimension d; we subsequently specialize them 
to the subcase (d = 1 and) d = 2. 

Let us make a comparison with the existing literature on the subject. The total 
energy and the forces on the sides of a rectangular box have been discussed in a 
lot of works, some of them using zeta regularization; here we only mention some 
of them. The foremost computation was performed by Lukosz [20; 2IJ for the elec¬ 
tromagnetic field, by means of exponential regularization and Abel-Plana formula; 
the same technique was used by Mamaev and Trunov [22, [23J (also see [2H 25]) to 
discuss, amongst other models, the case of a conformal scalar field 0. Alternative 
derivations of the total energy for a scalar field based on global zeta regularization 
were given by Ruggerio, Vilanni and Zimerman [an S2] in two and three spatial 
dimensions, and by Ambjprn and Wolfram [3] (see also [3j for the electromagnetic 
case) in the case of a multidimensional rectangular cavity for several boundary con¬ 
ditions. The same configurations were later re-examined by X. Li, Cheng, J. Li and 
Zhai [19] by means of a zeta strategy, and by A. Edery [HI E] using a so called “mul¬ 
tidimensional cut-off technique”. Let us also cite the papers by Estrada, Fulling et 
al. [15] OS], on which we return later in this Introduction. Finally, let us mention 
the monographies of Elizalde et al. unmu and Bordag et al. [5]; these can be taken 
as standard references for the study of global aspects. 

In the present paper we consider global observables, such as the total energy and 
forces, mainly to complete the analysis of local aspects (i.e., the VEV of the stress- 
energy tensor), which occupy most of our analysis. Our series representations for 
these global observables are different, but equivalent to the ones of [5] (which appear 
to converge exponentially, like ours). 

Concerning local aspects in the previous literature for a scalar field in a rectangu¬ 
lar box let us first mention the two seminal papers [1] and [3 by Actor. In these 

2 These authors also consider the cases of electromagnetic, (massless) gluon and spinor fields in 
(3 + l)-dimensional Minkowski spacetime; furthermore, they derive the renormalized average over 
the spatial domain of the stress-energy VEV. 
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papers d = 3, the framework is Euclidean and the author renormalizes by analytic 
continuation the effective Lagrangian density and the VEV (O|0 2 (a:)|O). The infor¬ 
mation contained in these works is equivalent, in the language of our papers, to the 
specification of the Dirichlet function D s (x, y) := (<5 X |(—A) s <5 y ) along the diagonal 
y = x (with A the Laplacian and x, y E fi), including its analytic continuation with 
respect to the complex parameter s. However, the VEV of the stress-energy tensor 
depends on D s (x, x) as well as on the derivatives of D s (x, y) with respect to x and 
y, at points of the diagonal (see Part I of the present series). For this reason the 
results of DEI are not sufficient to determine the stress-energy VEV, which in fact 
is not mentioned therein. 

The renormalized VEV of the stress-energy tensor is derived in a work of Svaiter et 
al. [30]; the latter employs, again, analytic continuation methods in a formulation 
closely related to the approach of Actor. However, in [30] the authors consider an 
infinite rectangular waveguide, rather than a box; more precisely, it is assumed d = 3 
and the spatial domain is H = (0,ai) x (0, a 2 ) x R Apart from the domain, there 
are other differences bewteen the approach of Svaiter et al. and ours; in particular, 
the methods employed in [30] ultimately yield a representation of the stress-energy 
VEV via series converging with polynomial speed. Even though this approach could 
be extended to treat a box domain (H = (0,ai) x (0,a 2 ) x (0,a 3 )), it would most 
likely yield similar series representations, whose polynomial convergence would be 
slower compared to the exponential convergence of ours. 

To go on, let us return to the already mentioned works by Estrada, Fulling et al. 
[T5lfl6] . Therein the configuration of a 2-dimensional rectangular box is analysed for 
both Dirichlet and Neumann boundary conditions (along with some related models, 
such as the Casimir piston and the so-called “Casimir pistol”). The cited works 
introduce a regularized version of the stress-energy VEV, based on an exponential 
cutoff; this can be expressed in terms of the so-called cylinder kernel, which was 
also considered for other reasons in Parts I and II of this series. The regularized 
stress-energy VEV is computed by the method of images, allowing to express it as 
a sum over infinitely many optical paths. By integration of the above VEV, the 
authors of [i5i rro] also obtain the regularized total energy and the regularized force 
acting on a side of the boundary. Their position of principle is that the theory with 
a cutoff is a more realistic description of the physical system under investigation; 
nonetheless, they also point out that the VEV of the observables for the system under 
analysis can be renormalized retaining only their regular parts with respect to the 
cut-off (an idea somehow related to what we call the “extended zeta approach”). 
Renormalization along these lines is carried out for global observables like the total 
energy, and hinted at for local observables (namely, for the energy density). 

The present paper is organized as follows. In Section [2] as in the other papers of 
this series, we report a summary of results from Part I to be used in the present 
work. Section [3] and the related Appendices |A] [B| |C] are the core the paper; therein 
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we treat the box configuration 0 = (0, eti) x ... x (0, ad) in arbitrary spatial dimension 
d. for Dirichlet boundary conditions. Our starting point is the heat kernel (<5 x |e _t4 <5y) 
for which we derive two different series representations capturing, respectively, the 
behavior for small and large t. As emphasized in Part I, the heat kernel determines 
an integral representation of the Dirichlet function ZT,(x, y) := (<5 X |(—A) s 5 y ), which 
can be used to construct the analytic continuation in s of the latter. Combining 
these general facts with the series representation for the heat kernel cited above, we 
ultimately produce series expansions for the analytic continuations of Z) s (x, y) and 
its derivatives, with the previously mentioned exponential speed of convergence; 
these determine all the local or global renormalized observables indicated before, 
from the stress-energy VEV to the force on each side of the box. 

In Section |4j the previous general results are specialized to the cases with d G {1,2}. 
For d — 1, we recover from a different viewpoint the results already obtained in 
Section 6 of Part I where, as a first example of our general formalism, we discussed 
a massless scalar field on a segment (0, a) (see this section of I for some references on 
this case, especially m ). For d = 2, working on = (0, cq) x (0, 02 ) and using the 
previously mentioned series expansions, we produce several graphs; some of them 
refer to the components of the stress-energy VEV for some choices of ai,a 2 , while 
the others are about the further observables in which we are interested. 

Before moving on, let us mention that many of the results presented in this paper 
have been derived with the aid of the software Mathematica for both symbolic and 
numerical computations. 

2 A summary of results from Part I 

2.1 General setting. Throughout the paper we use natural units, so that 


( 2 . 1 ) 


c = 1 , h — 1 . 


Our approach works in (d + l)-dimensional Minkowski spacetime, which is identified 
with R d+1 using a set of inertial coordinates 



( 2 . 2 ) 


X = 


the Minkowski metric is ( 77 ^) = diag(—1,1,..., 1). We fix a spatial domain O C R c/ 
and a background static potential V: —> R. We consider a quantized neutral, scalar 

field (p : R x 0 —> £ sa ($) (S is the Fock space and £ sa (S) are the selfadjoint operators 
on it); suitable boundary conditions are prescribed on dfl. The field equation reads 


0 = (~d a + A - V (x))0(x, t) 


(2.3) 


(A := y~)f_| da is the d-dimensional Lapla.cian). We put 


A:=-A + V , 


(2.4) 
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keeping into account the boundary conditions on dfl, and consider the Hilbert space 
L 2 {Q) with inner product (f\g) f Q dx f(x)g(x). We assume A to be selfadjoint 
in L 2 (fi) and strictly positive (i.e., with spectrum cr(*4) c [e 2 , +oo) for some £ > 0). 
We often refer to a complete orthonormal set (F k ) kefC of (proper or improper) eigen¬ 
functions of A with eigenvalues (u%) ke jc (u> k ^£ for all fce/C). Thus 

F k : —> C; AF k = c o k F k ; 

(F k \Fh) = 5(k, h ) for all k, h G K, . 

The labels k G K, can include both discrete and continuous parameters; J ' K dk indi¬ 
cates summation over all labels and 5(h,k) is tl ie Dirac delta function on /C. 

We expand the held 0 in terms of destruction and creation operators corresponding 
to the above eigenfunctions, and assume the canonical commutation relations; |0) G 
5 is the vacuum state and VEV stands for “vacuum expectation value”. 

The quantized stress-energy tensor reads (£ G R is a parameter) 

T I1U := (1 - 2 £) <9^0 0 ^ 0 - Q - 2 ^ ri^d^dx^ + VcJ?) - o 3^$) ( 2 . 6 ) 

in the above we put A o B := (1/2 )(AB + BA) for all A, B G C sa ($), and all the 
bilinear terms in the held are evaluated on the diagonal (e.g., d o d u (j) indicates 
the map x i-G d^cp^x) o d v (p(x)). The VEV (0|7j t! ,|0) is typically divergent. 

2.2 Zeta regularization. The zeta-regularized field operator is 


:= (K- 2 Ay u/i $ , 


(2.7) 


where A is the operator (2.4), u G C and k > 0 is a “mass scale” parameter; note 


that 4> u \ u =o — 0, at least formally. The zeta regularized stress-energy tensor is 


T; v := (l-2O^0 u o^0“-^--2eJ^^>^ A 0“ + V(0“) 2 j -2e0“o^0V (2.8) 

The VEV (0|T“ |0) is well defined for large enough (see the forthcoming subsec¬ 


tion 2.5); moreover, in the region of definition it is an analytic function of u. The 


same can be said of many related observables (including global objects, such as the 
total energy VEV). 

For any one of these observables, let us denote with F{u) its zeta-regularized version 
and assume this to be analytic for u in a suitable domain Uq C C. The zeta approach 
to renormalization can be formulated in two versions. 

i) Restricted version. Assume the map U 0 —)• C, F{u) to admit an analytic 
continuation (indicated with the same notation) to an open subset U C C with 
U 9 0 ; then we define the renormalized observables as 


F ren ■ *^”( 0 ) . 


(2.9) 
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ii) Extended version. Assume that there exists an open subset U C C with Uq cU, 
such that 0 GW and the map uEUq 4 T'(m) has an analytic continuation to 7/\{0} 
(still denoted with F). Starting from the Laurent expansion F(u) = J2k=-oc -^kU k , 
we introduce the regular part (RPF)(u) := Y^k=o^kU k and define 


Fren '■= (RP F)( 0) . 


( 2 . 10 ) 


Of course, if F is regular at u = 0 the defnitions (2.9) (2.10) coincide. 

Differently from the other papers of this series, the observables considered in the 
present Part IV will never exhibit singularities at u = 0; thus, we will never need to 
use the prescription (ii). 

In the sequel we will often refer to the stress-energy VEV, for which the prescription 
(i) gives 


(0\T^(x)\0)ren ■= (01T“„ (x) 10) 


u=0 


( 2 . 11 ) 


2.3 Conformal and non-conformal parts of the stress-energy VEV. These 
are indicated by the superscripts (0) and (■), respectively; they are defined by 


(0\f Ml/ \0) ren = (0\T^\0) ren + (£-&) (0\T^\0) ren , 

where we are considering for the parameter £ the critical value 


£d ■ = 


d-1 

~Jd 


( 2 . 12 ) 


(2.13) 


2.4 Integral kernels. If B is a linear operator in L 2 (f2), its integral kernel is the 
(generalized) function (x,y) e flxO 4 £>(x,y) := (S x \BS y ) (<5 X is the Dirac delta 
at x). The trace of B, assuming it exists, fulfills Tr B = f n dxB(x,x ). 

In the following subsections A is a strictly positive selfadjoint operator in L 2 (f2), with 
a complete orthonormal set of eigenfunctions as in Eq. (2.5). In typical applications, 
A is the operator (2.4). 


2.5 The Dirichlet kernel and its relations with the stress-energy VEV. 

For (suitable) s G C, the s-th Dirichlet kernel of A is 


r au _ 

D s (x,y) :=A~ s (x,y) = / F k (x)F k ( y) . (2.14) 

Jk u k 


If D is a bounded domain with smooth boundary dfl and A = —A + V (with V a 
smooth potential) is a strictly positive operator on L 2 (f2), considering the closure 
Q U dfl one can show that the map D s ( , ):flxfl4C, (x,y)4 D s (x,y) 
is continuous along with all its partial derivatives up to order j G N, for all s G C 
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with Jhs > d/2+j/2 ; moreover, the eigenfunction expansion in Eq. (2.14) converges 
absolutely and uniformly on with all the derivatives up order j for > d + j/ 2. 


Recalling Eq. (|2.8|), the regularized stress-energy VEV can be expressed as follows: 

. (2.15) 


<0|I£(x)|0) = «’ 


7+8 B»=i( x ,y)+ t-{ (8'8j+l'(i))Dii(xj) 


y=x 


= K 


(o|f”(x)|o) = (o|y"(x)|o) = o, 

(0|2)“(x)|0) = <0|7”(x)|0> = 

j -?)5ij(D»=i( x ,y) - (a I 'a sl +12(x))D^i(x,y)0 + 


(2.16) 


(2.17) 


2 C) f &x i xi ) ^\ 1 (^> y) 


y=x 


((01 T^ v (x)10) is short for (0|T“„(f,x)|0); indeed, the VEV does not depend on t). 
Of course, the map 0 —* C, x H> (0|T“ !/ (x)|0) possesses the same regularity as the 
functions x G 0 i —> Du±i (x. x). d^„Du +1 (x. x) (z,w any two spatial variables); so, 

due to the previously mentioned results, x H > (0|T^(x)|0) is continuous on 0 for 
> d + 1. 

Assume that Du=i and d zw Du+i (for z,w any two spatial variables) have analytic 
continuations regular at u = 0 (which happens for the configuration considered in 
this paper); then, we can define 


£> ± i(x,y) := D«±i(x,y) 


It—0 


d zw D i(x,y) := d zw Du+ i(x,y) 


tA=0 


(2.18) 

(2.19) 


The renormalized stress-energy VEV (0|T AtI/ (x)|0) ren := (0|T“ !/ (x)|0)| u=0 (see Eq. 
(2.11)) can be expressed as 


(0|Too(x)|0) ren = 


^W(x,y)+Q-£) {dfdyt+V (x))D + 1 (x, y) 


y=x 


(0|f 0i (x)|0) ren = (0|f j0 (x)|0) ren = 0 , 
(01 Tij (x) 10) ren = (01 Tji (x) 10) ren = 
i(x,y) - (<9 x ^ + V(x))L> + i(x,y)) + 

1 


, ( 2 . 20 ) 
( 2 . 21 ) 


( 2 . 22 ) 


+ ((2 -£) d x*yi -id x i x AD + i(x,y) 


y=x 


(More generally, in presence of a singularity at u = 0 we should consider the regular 
parts of the functions k u Du- 1 , K u d zw D „.+\). 















2.6 The heat kernel. For t G [0,+oo), this is given by 


K{ t;x,y) :=e u (x,y) = [ dk e tuj * F k (x)F k ( y) . (2.23) 

JK, 


If f2 is a bounded with smooth boundary dVt and A = — A + H (V smooth) is strictly 
positive, the map K( t; , ) : fl x fl —* C, (x, y) i-> Lf(t;x, y) is continuous along 
with all its partial derivatives of any order for all t > 0; moreover, the eigenfunction 
expansion in Eq. (2.23) converges absolutely and uniformly (the same holds for all 
the corresponding derivatives). 


2.7 The Dirichlet kernel as Mellin transform of the heat kernel. For suit¬ 
able values of sGC (see Part I), there holds 


A,(x, y) 


r»+C50 


F(s) 


di t s_1 iF(t; x, y) . 


(2.24) 


2.8 The case of product domains. Factorization of the heat kernel. Let 

A := — A + V and consider the case where 

O = hi x 0 2 9 X = (xi,x 2 ) ,y = (yi,y 2 ) , (2.25) 

P(x) = E 1 (x 1 ) + E 2 (x 2 ) (2.26) 

(fla C R da is an open subset, for a G {1,2}; di + d -2 = d): assume the boundary 
conditions on to arise from suitable boundary conditions prescribed separately 
on and <9f2 2 so that, for a = 1,2, the operators 

Aa ■= -A„ + V(x 0 ) (2.27) 

(with A a the Laplacian on Q a ) are selfadjoint and strictly positive in L 2 (f2 a ). Then, 
the Hilbert space L 2 (Q) and the operator A can be represented as 

L 2 {VL) = L 2 (Hi) <8> L 2 {VL 2 ) , ^ = + (2.28) 

This implies, amongst else, that the heat kernels K( t; x, y) := e tA (x, y), K a ( t; x a , y a ) 
:= e t ' Aa (x a , y a ) (a = 1,2) are related by 

JF(t;x,y) = A'i(t;x 1 ,y 1 )/F 2 (t;xi,x 2 ) . (2.29) 

Similarly, writing K( t), K a ( t) (a G {1,2}) for the heat traces of A and A a (a G 
{1,2}), respectively, we have 

K(t) = Ki(t) K 2 (t) . (2.30) 

The above considerations have obvious extensions to product domains with more 
than two factors. 
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2.9 Pressure on the boundary. This is the force per unit area produced by the 
quantized held inside 12 at a point x G <9f1 We first consider, for 3 ftu large, the 
regularized pressure p“(x) with components 

p”(x) := (0|T”(x)|0K(x) ; (2.31) 

here and in the remainder of this paper, n(x) = (n*(x)) is the unit outlier normal 
at x G <912. For Dirichlet boundary conditions the above definition implies 


PiW = 


7 Sij d x 'd y e + ^ d x i yj ) D u±i (x, y) n 3 (x) 


y=x 


We can define the renormalized pressure by analytic continuation as 


pr(x) :=P?(x) 


u =0 


(2.32) 


(2.33) 


(more generally, one should take the regular part if u = 0 is a singular point; this 
will never happen in the present paper). Alternatively, we could put 


K en (x):= , lim <0|T^(x')|0> r 

x'ES2,x —>-x 


n 3 (x) 


(2.34) 


Prescriptions ( |2.33[ ) ( |2.34[ ) do not always agree (for a counterexample, see Section 
5 of Part If 


of this series of papers). In Part I we conjectured that the two 
approaches agree when both of them give a finite result; this conjecture is true in 
the present case of a box as shown in subsection 3.6). 


Contrary to the previous works of this series [121 IB] [2] , here we first discuss the 
pressure and next pass to the total energy; we make this choice because quite dif¬ 
ferent computation techniques are employed for the stress-energy VEV and for the 
pressure, on the one hand, and for the total energy and integrated force, on the 
other hand. 


2.10 The total energy. The zeta-regularized total energy is 


S u := [ dx (0|T“ (x)|0) = E u + B u ; (2.35) 

Jn 

the second equality is proved after defining the regularized bulk and boundary en¬ 
ergies, which are 


K 


K 


E u := — / dx Du~i(x, x) = — Tr A 2 , 


B u :—k u [-—£) [ da(x) —' — fl»+i (x. v) 


'on 


dn y 2 


y=x 


(2.36) 

(2.37) 
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One has B u = 0 for 0 bounded and either Dirichlet or Neummann boundary con¬ 
ditions on <912. 


Assuming the functions (2.36) (2.37) to be finite and analytic for suitable u G C, 


we define the renormalized energies by the restricted (or generalized) zeta approach; 
for example, we put 

E ren := E u (2.38) 


1Z=0 


(again, if a singularity appears at u = 0 we should take the regular part). 

2.11 Forces on the boundary. Let £) C <912; we consider (for 3ftw large) the 
regularized integrated force acting on £), i.e., 


do : = / da(x) p“(x) 

Jo 


(2.39) 


(p“(x) is the regularized pressure of Eq. (2.31)). We can define the renormalized 
total force on £> as 


rgren ._ ^ 

on •— O; 


D 


11=0 


(2.40) 


(once more, the regular part should be taken if u = 0 is a singular point, a variation 
which we shall never employ in the present paper). Alternatively, we could put 


d r 0 en := / da(x) p ren (x) , 


(2.41) 


<o 


(pren(x) j s the renormalized pressure, dehned according to either Eq. (2.33) or Eq. 


(2.34)). In general the alternatives (2.40) (2.41) give different results. 


3 The case of a massless field inside a box 

3.1 Introducing the problem. In the present section we analyse the model of a 
massless scalar field confined within a d-dimensional box, with no external potential; 
more precisely, we assume 

12 = xf =1 (0, a,j) with cq > 0 for i G {1,..., d} , V — 0 . (3.1) 

The boundary <912 of the spatial domain is composed by the sides 

7 r p ,A := {x G R d I x p = A a p , x 1 e [0, a*] for i ^ p} 
for p e {1,..., d}, A e {0,1} ; 

for the sake of simplicity, we restrict attention to the case where the field fulfills 
Dirichlet boundary conditions on each one of these sides, meaning that 

<p(t, x) = 0 for all t G R, x G 7r P)A (p G (1,..., d}, A G {0,1}) . (3.3) 
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As a matter of fact, all the results to be reported in the following could be generalized 
to the case of Neumann or periodic boundary conditions, possibly including cases 
where different boundary conditions are prescribed on different sides of the box; 
moreover, the methods to be presented could be adapted with little effort to deal 
with the cases of a massive scalar held (V = m 2 ) and of a slab configuration (see 
Part I) where 

n = Qi x R d2 , fil = xti(0, Oi) c R dl (d 1 + d 2 = d) . (3.4) 


None of these generalizations will be considered in the present paper. 

3.2 The heat kernel. Similarly to the model with a harmonic potential consid¬ 
ered in our previous work [13] , in the present setting we are dealing with a product 


domain configuration, in the sense of subsection 2A Working in standard Cartesian 
coordinates (:r*)i=i,...,d, the Hilbert space and the fundamental operator A := —A 
can berepresented as 


A — A\ 
A; '■ = 


L 2 (n) = ($L 2 ( 0,Oi) , 

i= 1 

l <g>... <g> l +... +1 <g>... i 

d 2 


dx 2 


on L 2 (0, aj) 


... ® 1 ( 8 ) A d , 

11 e {!,...,<(}) i 


(3.5) 


(3.6) 


for each operator Ai we assume the induced Dirichlet boundary conditions in x l = 0 
and x z = at. 

According to the general considerations of subsection A8] (see, in particular, Eq. 
(2.29)), in this situation the heat kernel associated to A factorizes; more precisely, 
we have 


6'('t: x. y j - J1 Wt: x\ y 1 ) , 


(3.7) 


i =1 


where, for i G {1, Iij indicates the heat kernel of Ai . Hereafter we compute 

the heat kernels K ,, giving for all of them two distinct representations; these are 
suited to describe, respectively, the behaviour of the kernels K{ for small and large t 
(in a sense to be made more precise in the following) (n). Clearly, each one of these 


representations yields, in turn, alternative expressions for the total heat kernel (3.7). 

3.2.1 First representation for the heat kernel (useful for large t). Let us 

begin noting that, for any i G {1,..., d}, a complete orthonormal set of eigenfunctions 
(Fk^ki&lCi fc> r A in L 2 (0,aj), with corresponding eigenvalues (o^.)*,.^, is given by 


F ki (F) := \ - sin fax 1 ), uj 2 k := k 2 
w ai 


for ki G /Cj = <J - 

CLi 


UiTC 


rti = 1,2,3,... 


(3.8) 


3 For further information about approximate evaluations of the heat kernel for a box configura¬ 
tion, see, e.g., [3 H3 HU- 
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Using the eigenfunction expansion (|2.23|), we obtain for the 1-dimensional heat kernel 
K, the following expression ^ • 


+OO 


K. 


,, i is z i\ • P 2 ^ 

(t ;x\y) = —} e * sin - y ) sm I- x ). 

a.i V a,- / V a,- / 


(3.9) 


n,=l 


The above relation, along with Eq. (0, yields 

2 d 


K( t;x,y) = - V e-""* C„(x,y) , 

ai...a d ^ 


neN d 


where, for the sake of brevity, we put 

N:= {1,2,3,...} , 

2_2 d 

4==E" 


(3.10) 


(3.11) 


UJ~ ■■= > T ^ r , C n (x,y) J^sinT—y l )sin( — a; 1 ) for n := . 

i=i a * i=i k ai 


CI 7 , 


Eq. ( |3.10[ ) is easily seen to give a large t expansion for the heat kernel K (t; x, y) of 
A ; with this, we mean that the series over n G N c/ in the cited equation is mainly 
determined by the terms corresponding to small values of rq, for i G {1,..., d}. 

3.2.2 Second representation for the heat kernel (useful for small t). Let 

us move on and note that, with little effort (namely, writing the sines in terms of 


complex exponentials), we can rephrase Eq. (3.9) as 




a ■ 

e ' 


JxAix'-y*) J-£r- (x'+y*) 
e t — e * 


m=—00 


Using the Poisson summation formula (FI) and noting that 


'* + 00 7T 2 . 2 „ 

, i — ((x l ±y l )-2a.ihi)z a i - 

dz e a i e “i u y> ' =—= e 4 ‘ 


for /ijG Z with Z := (0, ±1, ±2,...} , 


Let us mention that the series in Eq. (3.9) could be explicitly evaluated to yield 




e~ a h- e 4t oJ-i 

2t ’ ) VTT 3 V 


(3.12) 


(3.13) 


. a t {x l +y l ) _°1 


2t 


where 6 * 3 ( , ) denotes one of the Jacobi elliptic theta functions; see [ 6 j for more details. 

5 The Poisson summation formula states that, for any sufficiently regular function / : R —>• C, 
there holds 

+00 +00 <-+00 

E = E /W > where />) ; = / e~ 2inhz f(z) 

n =— 00 h——oo ~ 00 


(/ is, essentially, the Fourier transform of /). Eq. (3.13) gives / in the case we are considering. 
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it follows from Eq. (3.12) that 

1 +00 


K l {i-x\y i ) = 


a/ 47 Tt 


£ 


( 2 aih i -(x % +y t )Y 
4t 


(3.14) 


hi ——oo 

The above identity can be rephrased as follows: 

1 




(3.15) 


where 


4 := 


1 for li — 1 

-1 for L — 2 ’ 


hi£ Z, /jG{l,2} 


H (I . 9 n , = , 17 *» '< = 1 

U, ‘ {X ’ V >- 1 fori, = 2 


(3.16) 


representation 


Eq. (|3.15|), along with Eq. (|3.7|), allows us to infer for the heat kernel the alternative 

(3.17) 


#(t;x,y) = 


(47rt) rf / 2 


5 > e " 


-f fe hi(x,y) 


hez d ,le{i, 2 } d 


where, for simplicity of notation, we have put (recall Eq. (3.16)) 

d 

h • i 1 • (li)i =5 ^1 • $li > 

d 

6 hl (x,y) := ^^{hi-UiXx 1 ,y l )f 


i =1 


(3.18) 


2 — 1 


Notice that, for small t, the sum of the series appearing in Eq. (3.17) is mainly 
determined by the terms corresponding to small values of |h,;|, for i G { 1 ,..., d}; thus, 
the mentioned equation yields a small t expansion for the heat kernel Ji(t;x, y). 

3.2.3 Considerations on the sign of some coefficients. Before moving on, 
let us emphasize a number of facts on the expansions (3.10) (3.17); we will resort 
to them in the following subsections, when performing the analytic continuation of 
the Dirichlet kernel and of its derivatives. 


b The very same result of Eq. (3.14) could be obtained via the method of reflections, starting 
with the heat kernel associated to the laplacian — on R, i.e. 


K(i\x,y) = 


1 


(47ft) d / 2 
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i) On the one hand, we have 


oj „ > 0 for all n G N" 


(3.19) 


ii) On the other hand, notice that 

&hi(x, y) ^ 0 for all h G Z d , 1 G {1, 2} d , x,y Gil 
In particular, since 

TT.Ui f [~2’H for/ * = 1 

[0,1] for li = 2 


(3.20) 

(3.21) 


(see the definition of Ui t in Eq. (3.16)), it follows that 

M x ,y) = ° 

( hi = 0, k = 1 , y l = x l G [0, a r ] 

for each i G {1, ..., d } , one has < or hi = 0, k = 2, y l = x l = 0 

[ or hi = 1, k = 2, y l = x l = di 

let us stress that this implies, in particular, 6hi(x, y) ^ 0 for x ^ y. 


(3.22) 


3.3 The Dirichlet kernel. Consider the integral representation (2.24) for the 
Dirichlet kernel D s in terms of the heat kernel K of A . We are going to construct 
the analytic continuation of D s in the style of Minakshisundaram (see BII)i to this 
purpose, let us fix arbitrarily 

Tg (0,+oo) (3.23) 

and re-express the cited integral representation as 


D.ix.vi = £>M (x, y) + £>)<> (x, y) , where 


£>< >) (x,y) : = 


q+oo 


r(s) 


dtt s K( t;x,y) , 


L)( <} (x,y) f di I s 1 K( t;x,y) 

r (s) Jo 


(3.24) 

(3.25) 

(3.26) 


(notice that D ( s >} and D'A’ depend on T, but their sum D s does not!). The idea 
we are going to pursue in the following is to substitute into Eq.s (3.25) (3.26), 
respectively, the large and small t expansions (3.10) (3.17) for the heat kernel K ([(). 


>«) 


' In these manipulations (and in some related computations) we often take for granted that 
certain series can be integrated or differentiated term by term. In all cases under analysis, rigorous 
justifications could be given using the Lebesgue dominated convergence theorem or the Fubini- 
Tonelli theorem, but we will not go into the details; the estimates on the series in Appendices |A| 
and [C| could be connected to such rigorous proofs. 
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3.3.1 Series expansion and analytic continuation of D^(x, y). Using Eq.s 
(13.101) (13.251), we readily infer 


Oi >) (x.y) = 


a 1 ...a d r(s) 


p +oo 

^ C n (x,y) / dt 

iSN d JT 


(3.27) 


Concerning the integral over (T, + oo), via the change of variable r := t (recall 
that > 0 for all n G N cZ ; see Eq. (3.19)), we obtain 

/•+oo /•+00 


dt t 8 - 1 e= u~ 2s 


olT 


dr T s - L e~ T = uj- 2s T(s,u z n T) forsGC, (3.28) 


where the last passage contains the upper incomplete gamma function (see 
p.174, Eq.8.2.2) 


“ + OO 


r(s, z):= dw e w w s 1 

J Z 

Summing up, we have 


(s G C, z G (0, + 00 )) . 


Df\^ y) = 


d\. ..&d r (s) 


W n 2 Sr (s^n T ) Cn(x,y) 


(3.29) 


(3.30) 


neN d 


The above expression can be used to evaluate derivatives of any order of the function 
Di > \x, y); for example, for any pair of spatial variables z,w, we obtain 

d zw D { >\x, y) = -- ^73 ^ u n 2sT ( s ^i T ) d zw C n (x, y) . (3.31) 


ai • ••(id r(s) 


neN d 


Let us anticipate that the series in the right-hand sides of Eq.s (3.30) (3.31) converge 
for all s G C, even for y = x (see subsection 3.4 and Appendices |A[ |C| for more 
details); so, Eq.s (3.30) (3.31) yield automatically the analytic continuations of the 
maps s 1 —y Di > \x, y), d zw D^ (x, y) to the whole complex plane. 

3.3.2 Series expansion and analytic continuation of D^(x, y). Proceeding 
similarly to what we did for the function Ds(x, y), we can use Eq.s (3.17) (3.26) 
to deduce 

r T 


£>M(x.y) = 


( 47r )d/ 2 F( s ) 


£ 

rt 1.— , 1 nl ^ do 


dt t s ~i~ l e~-i b ^y) 


(3.32) 


hez d ,ie{i,2} d 

To go on, for /3 ^ 0 and s G C satisfying the conditions in the forthcoming Eq. 
(3.34), let us introduce the function 


V s {/3) : = / drr 8 x e - ; 


(3.33) 
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it is easy to check that 


W) = 


s 1 for (3 = 0, > 0 


(3.34) 


/3 s r(-s,/3) for /? > 0 

where, again, T( , ) denotes the upper incomplete gamma function. Starting from 
Eq. (3.33), we find as well that 

d%P a {P) = {-l) l V s -z(P) for all i £ { 0 , 1 , 2 ,...} . 


(3.35) 



^ <} (x,y) = 


(4vr) d / 2 r(s) 


£ * ?.-* 


hez d ,ie{i, 2 } d 


fchi(x,y) 

T 


(3.36) 


The above result can be used, along with Eq. (3.35), to infer analogous representa¬ 
tions for the derivatives of any order of for example, if z, w are any two spatial 
variables, derivating term by term Eq. ( |3.36| ) we obtain 

d zw D ( <\^,y) = 


T s ~i~ 2 


(4 7r )d/ 2 E( s ) 


£ * 


hez d ,ie{i, 2 } d 


^s-i-2 


h ~Y )d z b h id w b hi 


TV s _d_ i 

b 2 i 


) d zw b hl 


(3.37) 

x,y) • 


Let us point out that, due to the results reported in Eq. (3.22), we have 

&hi( x , y) = 0 only for y = x and for 
a hnite number of terms in the series of Eq.s (3.36) (3.37). 


(3.38) 


The above mentioned terms of Eq.s (3.36) (3.37) deserve special attention, and we 
must use the Erst equality in ( |3.34[ ) to evaluate them; on the contrary, for the in¬ 
finitely many terms with &hi( x j y) > 0 , the second line in Eq. (3.34) gives expressions 
in terms of upper incomplete gamma functions. In this way we obtain 

/ \ 


D^(x,y) = 


(47r) d /2 T(s)(s — f) 


£ * 


hez d ,ie{i, 2 } 

\s.t. 6 h i(x,y)=0 / 


+ 


(3.39) 


(47r) d / 2 T(s) 


£ 


hezyieji^E 
s.t. fe h i(x,y)>0 



r(j-», y ]](*.y) 
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Let us repeat that the first sum in the above expression contains finitely many terms. 


Notice that the term in the first line of Eq. (3.39) is related to the first equality in 


Eq. (3.34) which, in principle, would require JLs > d/2; however, this term makes 
sense for all complex s except s = d/2 , where a simple pole appears. Moreover, the 


series in the second line of Eq.s (3.39) can be proved to converge for any complex 


s; we defer a comprehensive discussion of this statement to subsection 3.4 (see also 
Appendix [A]) . 


In view of the above remarks, Eq. (3.39) gives automatically the analytic continua¬ 


tion of Dg(x, y) to a meromorphic function of s on the whole complex plane, with 
a simple pole singularity only at 

s = d/2 (3.40) 


y = x due to Eq. (3.22). 


for x, y such that the first sum in Eq. (3.41) is non-empty; this happens only for 


A similar analysis can be made for the derivatives of D^. For example, if z,w are 


any two spatial variables, we obtain the following expression from Eq. (3.37): 


d zw D^\x,y) = 


(4n) d / 2 r(s)(s - | -1) ( £ /.^Mx,y)) + 

v ' wv 2 ’ \hez d ,ie{i,2} d 


(4 7r )d/2p( s j 


E * 


hez d ,ie{i, 2 } d 

s.t. b h i(x,y)>0 


s.t. 6 h l(x,y)=0 

b hi 2 2 frf^+2-5,^ ]d z b h id w b h i + (3-41) 


- r( d i+ls, b f )Kid zw b hl 


(x,y) 


Again, the first of the above two sums is made of finitely many terms; besides, 
contrary to what one conld expect from Eq. (3.37), this sum contains no term with 
the first order derivatives <9<„6hi( x ) y), <9^hi( x jy) because they vanish if &hi( x jy) = 
0. By considerations analogous to those made above for Eq. (3.39) (based on 
the convergence of the second sum for all s G C), we infer that Eq. (3.41) gives 
automatically the analytic continuation of d z 




(x, y) to a meromorphic f un ction 
of s on the whole complex plane, with a simple pole singularity only for y = x at 


s = d/ 2 + 1 


(3.42) 


Let us stress that Eq.s (3.39) (3.41) are just the original Eq.s (3.36) (3.37), rewritten 


separating the terms with &hi( x jy) = 0 for a better understanding of the behaviour 
with respect to s . In the sequel, even when considering meromorphic continuations, 


we will always refer to the more concise representations (3.36) (3.37). 
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3.3.3 Conclusions for the Dirichlet kernel. Using Eq. (3.24) and expressions 


(3.30) (3.36) for the functions D^, respectively, we obtain the analytic con¬ 
tinuation of the full Dirichlet kernel D s (x, y) to a meromorphic function on the 
whole complex plane. Similar results hold for the derivatives of the Dirichlet kernel 


see Eq.s (3.31) (3.37)) 


The only singularity of D s (x, y) is a simple pole for 

y = x and s = d/2 , 


(3.43) 


while d zw D s (/x.,y) (for any pair of spatial variables z,w) has a simple pole for 


y = x and s = d/2 + 1 . 


(3.44) 


In particular the analytic continuation of Dy^i (x, y)| y=x and d zw Du± i(x, y)| y=: 


y=x; 


required for the evaluation of the regularized stress-energy VEV and pressure, are 
both regular at u — 0 . 


3.4 Convergence and remainder estimates for the series in Eq.s (3.30) 


(3.31) and (3.36) (3.37). This subject is discussed in more detail in Appendix 
|A[ here we only report the main results. In the mentioned Appendix we show that 
the series cited in the title of this subsection are absolutely convergent; moreover, 
we derive fully quantitative remainder estimates when these series are approximated 
by finite sums. 

To report here these estimates we need some notations, introduced hereafter. First 
of all we put 


a := min {a*} , A := max {a*} ; 


z := 


d \V2 
„2 


£4 


(3.45) 


for z = n G N d or z = h 6 Z d . 


, i= 1 


Besides, for N E (2 y/d, +oo), aE (0,1), /3 E (0, +oo), a, p G R, we set 

HpPicx, (3]cr, p) := (3.46) 


7T 


d/2 


N-Vd 


,d -1 


T(a, (l-a)pN 2 ) T 


(1 -clY{clP)*¥t{$) \N-2Vd 

it can be shown that there holds the asymptotic expansion 
yj-d/2 2g—4a/3d 


( d ^]a(3{N-2Vd) 2 \ ; 


V 2 




a(l-a)T(f) 


e -f3N(N-4apVd)jy-2a+p+d-4 ^ _|_ Q^]\[~ 1 ^ 


(3.47) 


for N +cx) 
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Finally, we put 


c i%.(<r> N ) '■= max 




v 2cr / / \ \2cr 


(3.48) 


Having introduced the above notations, in the next two subsections we report the 


«) 


and 


remainder estimates of Appendix [a] for the series expansions of D^, D 
of their derivatives. In all cases the remainder is controlled by the function H 


due to the exponential decay of this function for large N (see Eq. (3.47)), good 
approximations of all the series under investigation can be obtained by just summing 
the first few terms. 


3.4.1 Estimates for the series (3.30) (3.3l|). Consider any s e C; keeping in 
mind Eq. (3.30), for any N G (0, +oo) let us write 


£)(>) | 


^S( x >y) : = 


x ,y) := 


^ >} ( x >y) = ^S( x >y) + ^ifi( x >y), 

2 d 


a 1 ...a d T(s) 


a 1 ...a d T(s) 


W n 2Sr ( S ^n T ) ^(x, y) , 


neNt |n|<iV 


(3.49) 


W n 2 ' r ( S ^n T ) ^nfoy) • 


iGN d , |n|>JV 


This equation implies a similar representation for the (analytically continued) deriva¬ 
tives of in terms of the derivatives of and R[^. For the remainder function 

i?^(x, y) and for its derivatives with respect to any two spatial variables z,w, we 
have the following uniform estimates: 


#S( x ,y) 




max a 


25R s ^2SRs\ 


ai...a (i 7r 25Rs |r(s)| 


h n ] f a, ^ ; 5 is, -2fts 


’ A 2 


for either Jhs^O, N >2y/d or U?s<0, N>2y/d+ — 

7r V a T 


(3.50) 


d zw R^(-x, y) 




max a 


2 JRs 


^ 2 JRs 


ai...a d a 2 7r 2 i SRs_1 )|r(s)| 


H‘i? I a, -Us, 2{l-» S ) 


’ A 2 


(3.51) 


for either 1, N>2\[d or !Rs<l, N>2y/d+ — \ — —4— 

7r V aT 

In the above, a is a parameter that can be freely chosen in (0,1); of course, the best 


choice is the one minimizing the right-hand sides of Eq.s (3.50) and (3.51), which 


depends on the other parameters (e.g., N,T ) involved in these considerations. 
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3.4.2 Estimates for the series (3.36) (3.37). Let s G C, and exclude the case 
(3.43); keeping in mind Eq. (3.36), for any N G (0, +oo) we put 

y) = L>^(x,y) + i?^(x,y) , 

£*S( x >y) : = 


^S( x >y) : = 


T s ~i 

E * v -i 1 

^6hi(x,y) 

(4vr) rf / 2 r(s) 

V T 

hez d ,|h|^v, ie{i, 2 } d 


T s ~-2 

E 

hez d , |h|>iv, ie{i, 2 } d 

^M x ,y) 

(47r) d / 2 r(s) 

v T 


(3.52) 


The above equation can be used to derive similar representations for the (analyt¬ 
ically continued) derivatives of D in terms of the derivatives of and 

with the exclusion of the case (3.44). For the remainder function and for its 
derivatives with respect to any two spatial variables z\ w J (■ i,j G {1,..., d}), we have 
the following uniform estimates: 


R ijU x ’ y) 


C^s-l,N) w ( d 


Vd\2 




7T d / 2 |r(s)| 


H7 a, 


;--3fo,2Ks-d 


d d 1 (^ts--)T 

for either $ls^-,N>2Vd or N>3Vd+-\ --— 

2 ' 2 a \ a 


(3.53) 


d z i w i y) 




(3.54) 


7T rf / 2 |r(s) I 


1 + Aj .4 2 C« (» S - d - -2. n)H^( a, 


- ;^+2-Ss, 2Sts-d-2 


“ ^ K ’ > ;- + l-Xs,2Xs-d-2\ 

T 2 


for either —+ 1, N> 2\fd or $ls>- + l,N>3Vd +-\^—— --— . 

2 2 a y a 

Again, the parameter a can be freely taken in (0,1) and it is conventient to choose 


for it the value minimizing the right-hand sides of Eq.s (3.53) (3.54), keeping into 


account the choices made for the other parameters (in particular, N). 

3.5 The stress-energy tensor. Consider the representations deduced in subsec¬ 


tion 3.3 for the analytic continuations of the Dirichlet kernel and of its derivatives. 
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Resorting to Eq.s (2.15 2.17), we obtain the following expressions for the components 


of the regularized VEV of the stress-energy tensor: 

<o|?;„(x)|o) = T£<»(x) + r;y>(x), 


(3.55) 


where, for • equal to > or <, T P y (x) has the espression corresponding to Eq.s 
(2.15 2.17), with D s replaced by D^\ Thus, for i,j 6 {1, we have 


->«,(•) _ 


L oo 


x — n 


J D% (x, y) + 0 - t\d x, d y iD { 'l (x, y) 


T“’ (#) (x) = T“’ (,) (x) = 0 , 


y=x 


■Oj 


L J'0 


T«' W (x) = K “ 


(j - d (pQ 1 (x, y) - (x, y)) + 


2 £ )dx l yl 


^d x i xj )D^ i(x,y) 


•u+1 1 
2 


(3.56) 

(3.57) 

(3.58) 


J y=x 


To proceed notice that, concerning the analyticity of the above functions, there hold 


considerations analogous to those presented in subsection 3.3.3 for the Dirichlet 


kernel and its derivatives; in particular, it appears that u = 0 is a regular point 
for each component of the regularized stress-energy VEV so that, according to the 
restricted version of the zeta approach, we can simply put 


(0|7V(x)|0}„„ := <0|T“ (x)|0) 


11=0 


(3.59) 


In the forthcoming subsections 4T and 4.2 , dealing with the cases d — 1 and d — 2, 
we will use approximate expressions for all the co mponents o f (0|T A1 ^(x)|0) r . en ob¬ 
tained replacing each Dirichlet function Dg in Eq.s (3.56 3.58) with the truncations 


of a fixed (sufficiently large) order N, given by Eq.s (3.49) (3.52). Let us re¬ 


call that we have explicit remainder bounds for these truncations 


see Eq.s (3.50) 


(3.51) and Eq.s (3.53) (3.54)); these will allow us to infer error estimates for the 
approximate expressions of (0|T Miy (x)|0) ren described above. 

3.6 The pressure on the boundary. We refer to the general discussion of sub¬ 


section 2.9 so, we have two alternative definitions for the renormalized pressure on 


each of the sides Tr Pt \ of the box (pG {1,..., d}, Ag{ 0,1}; see Eq. (3.2)). 

Let x be any point interior to one of the sides let us stress that we exclude x 
to be on an edge of the box (i.e., on the intersection of two or more sides), where 
the outer normal is ill-defined. As an example, let us assume x to be an inner point 
of the side 7T lj0 , so that the unit outer normal at x is n(x) = (—1,0, ...,0). We first 
consider the regularized pressure 


P“(x) := <0|7-(x)|0Mx) = - (0|T)“(x)|0) 


(3.60) 
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this can be expressed using the general rule (|2.32|) for the case of Dirichlet boundary 
conditions which, in the present case, gives 


K 


p“(x) = -5 n — d x i y iDu±i(x,y) 


y=x 


K 


- 5n 4 


d x iyiD^(x,y) +h>S(x,y) 


y=x 


y=x 


(3.61) 


The equality in the second line of the above equ ation , involving the derivatives of 
the Dirichlet functions D^h, follows from Eq. (3.24) and will be useful for later 

purposes. " - 

On the one hand, we can put 


r(x):=p“(x) =-(0|^(x)|0) 

u=0 


u =0 


8il 

T 


d x i y iD { >l(x,y) _ + d x i y i ( x , y) 


«)( 


y=x 


y=x 


(3.62) 


On the other hand, we have the alternative definition 


PT(x):=f , Jim (Olfi^xOlO)^^' 

\ x'eS2, x'—>x J 


= 


lim (0|Tji(x / )|0) ren ) . (3.63) 

x'eo,x'— 


In the next two paragraphs we show the equivalence of the alternative prescriptions 


(3.62) (3.63) and the non-integrable divergence near the edges of the renormalized 


pressure (evaluated, equivalently, according to either of the two cited prescriptions). 
Before proceeding to the proof of the above statements, let us anticipate that in 


subsections 4.1 and 4.2 (dealing with the cases d = 1 and d = 2, respectively) we 


will evaluate the pressure starting from Eq. (3.62) and substituting the functions 


D[% therein with the truncations D ( '*j 2 N of a sufficiently large order IV; the errors 
of these approximants will be evaluated using Eq.s (3.50) (3.51) and (3.53) (3.54). 


3.6.1 Equivalence of prescriptions (3.62) (3.63). The proof of this equiva¬ 


lence, given hereafter, uses arguments similar to the ones proposed in subsection 4.6 
of Part II, where an analogous statement was derived for a domain D bounded by 
orthogonal hyper-planes. 


8 In the application of Eq. (2.32) to the present case, we use the previous expression for n(x) 
and the fact that 

d x i y oD s fx.,y) = 0 for all i,j £ {1, ...,dj such that i ^ 1 or j ^ 1 ; 
y— x 

this follows straightforwardly from the Dirichlet conditions prescribed on the boundary of Q and 


from the eigenfunction expansion in Eq. (2.14|, taking into account the factorized structure of the 
eigenfunctions in our case. 


23 






























Let us consider again an inner po int x = (0, x 2 , x d ) of the side 7 Ti j0 . Whe n 
expressing (0|f^(x)|0 )| u=0 (see Eq. ( 3.62 )) or lim x /_> x (0|T i i(x , )|0) ren (see Eq. ( 3.63 )) 
in terms of the Dirichlet kernel and of its series expansions, there is only one type of 
potentially troublesome terms, which could give different contribut ions in the two 
cases. These are terms arising from the summand T((' (x) i n Eq. (3.58), when we 
use for D i < ' ) and its derivatives the series expansions (3.36) (3.37); more precisely, 
potential troubles could arise from terms in the cited expansions with 

&hi( x , x) = 0 and 6 hi( x/ ) x 0 7 0 for x 7 G , (3.64) 

corresponding to the choice 

hi — 0, li = 2 and hi — 0, h — 1 for i 7 1 . (3.65) 


Consider the expression for T i “’ (< ' ) (x / ) obtained using Eq. (3.58) and Eq.s (3.36) 
(3.37), where x' is either a point in the interior of the spatial domain or the 
boundary point x under consideration; the contribution to this expression from 
problematic terms, with (h, 1) as in Eq. (3.65), is 


fl \ k u T u 2 1 

u — d— 1 

\4 7 (47r) d / 2 r(^±4) 

- Vu-d-i (z) + zVu-d-i(z) 

2 2 2 


;1 =f( u , x>1 ) • (3-66) 


(x' 1 ) 

Z - 7jT~ 


To show the equivalence between prescriptions (3.62) (3.63) at the point x in con¬ 
sideration, we must show that 


f(u, 0 ) 


= lim 

u =0 WsfOjai), x' 1 - 


f(u,x 


/!> 


U = 0 


(3.67) 


where | u=0 indicates, as usual, the analytic continuation at u — 0. In order to check 
the equality in Eq. (3.67) we first notice that 


f(u, 0) = 6. 


ll 


\ k u T u 2 1 

l 

c 

p 

t-H 

1 

1 

3 

1 _ 

/ (4tt) d / 2 T(^) 

[2 2 J 


y, 


, , 1_A «" T ^ 

Jl 1 4 7 (47r) d / 2 r(^) 


(3.68) 


(the second equality relies on Eq. (3.34) for V s (0)). The last expression in (3.68) 
gives the analytic continuation of f(u, 0 ) to the whole complex plane; in particular, 
when evaluated at u — 0 , it yields 


f(u, 0) 


u =0 




(3.69) 
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On the other hand, returning to Eq. (3.66) we see that, for any x' 1 G (0, cq), 
f(u,x rl ) 

- z=(x' 1 ) 2 /T 


- (1 \ 

2 

d T1 _ , . _ , 


(4t rT)^ 

-— V d+i(z) + zT d+a(^) 

2 2 2 


(3.70) 


Resorting to the second equality in Eq. (3.34) for the functions V,(z) and using a re¬ 
cursive relation for the upper incomplete gamma function (see [29], p.178, Eq.8.8.2), 
we obtain 


f(u,x 


/!> 


= Sn T-e 

u= o \ 4 


uT) 2 


(4t r T)^ 


(3.71) 


comparing with Eq. (3.69), we immediately obtain the desired relation (3.67), that 
is f(u, 0)| u=0 = lim x ,i e(0;ai)ia ,,i^ 0 (f{u, . 

This concludes our analysis of the renormalized pressure at points in the interior of 
7Ti )0 ; needless to say, the very same results also hold for the pressure acting on any 
other of the sides delimiting 12 . 

3.6.2 Non-integrable behaviour of the renormalized pressure near the 
edges. Let us first remark that at points on the edges of the box (i.e., on the 
corners which appear whenever d > 1) the outer normal and, consequently, the 
pressure are both ill-defined. 

In this paragraph we discuss the behaviour of the pressure at points in the neighbor¬ 
hood of the edges; more precisely, we show that the renormalized pressure evaluated 
at inner points of one side diverges in a non-integrable manner when moving towards 
anyone of the edges. 

In order to fix our ideas, let us consider the renormalized pressure p^ en on the side 
7T 1)0 ; in the following we discuss the behaviour of this quantity in a neighborhood of 
the corner placed at x = 0, i.e., at the intersection of all the sides 7r Pi0 (p G {1,..., d}). 


To this purpose, let us consider the expression (3.62) for p \ en ; by considerations 


analogous to those of the previous paragraph, we see that contributio ns di verging 
for x —>■ 0 arise from d x i y iD^j} 2 when we use for it the series expansion (3.37). More 
precisely, these contributions arise from the terms for which 

frhi(0, 0) = 0 and 6hi(x, x) ^ 0 for x interior to 7Ti )0 , x^O ; (3.72) 

by a simple inspection, these terms are seen to correspond to the following choices 


hi = 0 for all i G {1,..., d} , 

U G {1, 2} for i G {!,..., d} and U = 2 for at least one i G {2,..., d} . 


(3.73) 


As an example, let us focus on one of the above terms; for any fixed p G {2, ...,<2}, 
we consider the one with 


hi — 0 for all i G {1,..., d} and l, — 


2 for i G {2, ...,p} 

1 for i G {p + 1,..., d} 


(3.74) 
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(<) 


With simple but long computations ( J ), the contribution of this term to d x i y i 
can be expressed as 


2 (-1 ) d ~ p 
(4vr T) 


d +1 ' — 


V 


d+1 


(. x 2 ) 2 + ... + (x p ) 2 


T 


= 9 P , i( x ) 


(3.75) 


Using the second relation in Eq. (3.34) and resorting to the asymptotic expansion 
(see J23], p.178, Eq.8.7.3) 


r(s, z) = r(s)(l + 0(z s )) for > 0, z —> 0 , 


we readily infer, for x —y 0, 

2 (-i)^r(^) 

0p,i( x ) =-TTTii+i- 

(47t)^ 


yd -\-1 


(1 + 0{z d+1 )) 


- Z=tJ(x 2 ) 2 +...+(xP ) 2 


(3.76) 


(3.77) 


The above expression shows the non-integrable divergence of g p j (for allpG {2,..., d}); 
the same conclusion holds as well for the terms corresponding to all the other choices 


of (h, 1) in Eq. (3.73). Let us stress that no cancellation of the divergent terms can 
occur, since all the contributions with the same degree of divergence happen to have 


the same sign; besides, due to the remainder estimate (3.54), no compensation of 


these terms can either arise from the full series expansion (3.37). 

The above comments prove the non-integrable behaviour of p[ e ”(x) for x —y 0; let 
us remark that this fact is of utmost importance when attempting to evaluate the 
total force acting on any side of the box, a topic we discuss in detail in subsection 

m 

3.7 The total energy. First of all, let us recall from subsection 


2.10 


that the total 


energy consists of the sum of both a bulk and a boundary contribution; since the 
latter vanishes identically due to the Dirichlet conditions assumed on the boundary 


(see the comments below Eq. (2.37)), we only have to discuss the bulk term. 
Consider the representation (2.36) of the bulk energy ( 10 ): 


K 


E u — / dx Du -1 (x. x) 


IQ 


9 In particular note that, with the present choices (3.73) (3.74), we have 

Mx,y) , „ = 


y=x, x 1 —0 


d x ib h i(x,y) 


y—x, x 1 =0 


= d y i6hi(x,y) 


= 0 


d x i v ib h i(x,y) 


y—x, x 1 —0 


y—x, x 1 —0 

10 Let us mention that, in order to evaluate the bulk energy E u , we could proceed in an alternative 
manner. Namely, we could consider the representation of E u in terms of the trace Tr A 2 " (see 


Eq. (2.36)) and determine the analytic continuation of the latter moving along the same lines we 
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Using the expression (3.24) for the Dirichlet kernel, we readily infer 

E u = E u ’ (> ) + U“’ (<) where 


2 


I dx 1 ... dx d -D 1 ^! (x, x) for •€{>,<}. 

'(0,a 1 )x...x(0,a d ) 2 


(3.78) 


In the next two paragraphs we discuss the analytic continuation of E u at u — 0, 
giving the renormalized bulk energy. To this purpose we introduce series expansions 
for the functions E u, ^ > ' > and E u, ^ < \ which can be used to build their analytic con¬ 
tinuations; we also give fully quantitative remainder estimates for these series (for 
more details, see Appendices [B| and [C|) . 

3.7.1 Series expansions and analytic continuations for E u, ( > \ E u, ^ < \ Con¬ 
sider the expression (3.78) for the regularized bulk energy E u ; hereafter we give series 
expansions for the two addenda E u ^ > ' ) and E u, ^ < \ ultimately yielding the analytic 
continuations of these functions to the wh ole co mplex plane. 

On the one hand, inserting the expansion (3.30) for into Eq. (3.78), we obtain 
(see Appendix [B| 


_£“>(>) — 


K 


2IW 


IT I] 


nEN d 


u- 


J U D T 


(3.79) 


The right-hand side of above equation can be proved to converge for all u G C 
(see the next paragraph and Appendix 0; thus, Eq. ( |3.79D gives the analytic 

continuation of E u ’( > > to the whole complex p lane, in particular at u — 0. _ 

On the other hand, inserting the expansion ( |3.36 ) for D i < ' ) into Eq. (3.78), after 
some effort we obtain (see again Appendix |BJ) 


£“>(<) — 


fiUTAO ^ (-l) d -P 


_ _V 

_2q A—/ (rj — ridnl 


A(J lP 




2 d + i r(^i) (d- p )\p\ ^ (ttT)p / 2 

y 2 ) p= o v yj f a&Sd \ > hgZP 


u—p—1 
2 


B «A h ) 

T 


(3.80) 


In the above S d indicates the symmetric group with d elements and we have put 
<1,0 := 1 , Z° := {0} , R,„fOi := 0 , 

(3.81) 


&(T,P * 


n 

i =1 


^cr(i) 


B ajP ( h) := y ^(aa(i)hj) 2 for aeS d , pe{ 1, ...,d} ; 


«=i 


note that the term with p — 0 in Eq. (3.80) is just (—l) d P™=i(0). 


followed for the Dirichlet kernel, using the heat trace /\(t) in place of the heat kernel /\(t;x, y). 
Nonetheless, since the small t expansion of the heat trace K{ t) for the present configuration involves 
quite cumbersome expressions, we prefer to avoid this approach. Another advantage of the methods 
proposed in this subsection is that, after minor variations, they also allow to evaluate the total 
force on the boundary (see subsection 3.8). 
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Let us stress that the functions V s in Eq. (3.80) must be evaluated according to Eq. 


(3.33); in particular, recall that the first relation in Eq. (3.34) gives 77,(0) = 1/s. 
Thus, for all p G {0, the terms in the series (3.80) with B ap (\\) = 0, i.e., 

those with h = 0 (see Eq. (3.81)), are singular at u = p + 1 where they have a 


simple pole. The series obtained from the right-hand side of Eq. (3.80) removing 


the finitely many terms with h = 0 is (rapidly) convergent for all u G C, a fact that 
we discuss in the next paragraph and in Appendix |C| 


Because of the above considerations, the expression (3.80) gives the analytic contin¬ 


uation of E u to a meromorphic function on the whole complex plane, with simple 
poles at 

u G {1,2,...,<2 + 1} • (3.82) 

Summing up, u = 0 is a regular point for the analytic continuations of both E u ^ > " > 
and E u, ( < > so that, according to the restricted zeta approach (see Eq. (2.38)), we 
can put 


E ren = E u ’ <y> " > 


u=0 


+ £“’ (<) 


u =0 


(3.83) 


where the two addenda on the right-hand side simply indicate the expressions (3.79) 


and (3.80) evaluated at u — 0. 


3.7.2 Estimates for the series (3.79) (3.80). The scheme followed in this 

paragraph closely resembles the one of subsection T4 where we dealt with the 
series expansions for the Dirichlet functions we retain here the same notations 


introduced therein (see, in particular, Eq.s (3.45) (3.46)). 


Let us consider the series in the right-hand side of Eq. (3.79), and the series in 


the right-hand side of Eq. ( |3.80| ) after removing from it the finitely many singular 
terms with h = 0. In Appendix [C] we prove the convergence of these series for all 
u G C, and derive the remainder estimates reported hereafter. For N G (0, +oo), 
let us write 


£“’(» = E u / >] + i^ (>) , 


J N 


K 


2V (—) ^ 

V 2 > n eN d ,|n|^iV 


R 


■“>(>)_ 

■N 


K 


2T(^1) ^ 

v 2 ' neN d ,|n|>jV 




'N > 



.ulT 

n V 2 


i-rf"' 1 

,uiT 

n \ 2 



(3.84) 


the remainder A/b has the bound (Appendix [c]) 


R 


N 


max(a s “ 1 ,A^ U ') t d \ ( n 2 T Ru — 1 
^ - . . J . Hy I a, — ; —— , l-Ku 


2 d + 1TT^- 1 |T(^)| 

for either ^ 1, N > 2 \fd or < 1, N > 2 Vd - \ - -^ . 

7r V 2 aT 


(3.85) 
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(3.86) 


Again for N G (0, +oo), let us put 

£“■«) = E^ <} + R u „ {<) , 

-i y- p v 


E 




N 


k’TV 

IX—1 ' 
2 


2 d + 1 r(^A 


E 


(d-py.pi. ^ 


a cr,p \ -r, 

— -g 2^ / «=!=! 

1 heZP, |h|^JV 


A 


«.(<) 


k“TV 


uE 


-i) 


d—p 


E 


*(T,P 


>r) 


E 


Vu-p-1 

2 


E a , p (h) 


Bg t p(h.) 

T 


2 "IW) ^ (d-p)!p! ^ v . hEZ , i|b|>Jf 

(where h| y A y + ... + / 1 ( j : the remainder /1’^- : which contains no singular term 

with h = 0, fulfills (see, again, Appendix [C| 


R 


“.(<) 

AT 






9fax 


2 d+ i |r(^i)| 


E 

p= i 


max(a s ” p 1 ,A 5R “ p x ) 
7 tp/ 2 (d—u)!n! 



or p+1 — 


, 57u—p— l) (3.87) 


for either 3fiw ^ 2, N > 2\fd or > d+1, N > 2'/d+- \ — --— . 

a V 2a 

As sketched in Appendix O we could give reminder estimates even for the case 2 < 
^ d + 1, excluded from (3.87); however, these would involve rather complicated 


expressions. Taking into account that, in the sequel, we will be mainly interested in 
the case u = 0, we prefer not to report these cumbersome expressions. 


In subsections 4T and T2, moving along the same lines as for the VEV of stress- 
energy tensor, we will use the truncations E and E of a fixed sufficiently 
large order N to obtain approximate expressions for the functions E u ^ and E u, ^ < \ 
respectively; these will be used to evaluate the renormalized bulk energy (see Eq. 


(3.78)), giving explicit errors estimates. 


3.8 The total force on a side of the box. Let us consider the framework of 


subsection 2.11; following the general scheme outlined therein for boundary forces, 
we can in principle consider two alternative approaches to define the total force 
acting on any side of the box. 

As an example, let us focus on the force acting on 7Ti i0 , he., the side contained in 
the hyperplane {x 1 = 0}; recall that the unit outer normal at points x interior to 
Tl,o is n(x) = (-1,0, ...,0). 

We first consider the regularized integrated force 


3T,o := / rf a(x) ™*( x )p“( x ) = ~ da{x) p“(x) 


(3.88) 


' 7T1,0 


' 7T1,0 
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(compare with Eq. (2.39), here employed with D = 7Ti j0 ), where p“(x) indicates the 


regularized pressure (3.60). We will prove in the sequel that this can be analytically 


continued up to u — 0, so that we can define the renormalized integrated force as 


r^ren . r*u 

Ol,0 Ol,C 


u =0 


(3.89) 


On the other hand, we have the alternative prescription (corresponding to Eq. 


(2.41), with D = 7r 10 ) 


:= / da(x) n ! (x)r(x) = - / da(x) p[ e 

J 7T1,0 J 7T1,0 


(3.90) 


where p^ en 
prescription (3.62) or (3.63). 


is the renormalized pressure, defined equivalently according to either 


As a matter of fact, we know from the previous subsection that the renormalized 
pressure diverges in a non-integrable manner near the edges of the box; in conse¬ 


quence of this, the prescription (3.90) gives an infinite value for the total force on 


7r li0 . Since this result is patently physically unacceptable, in the following we only 
consider the approach (3.88) (3.89). 


Let us therefore consider the regularized expression (3.88); using Eq. (3.61) for the 


regularized pressure on 7r 10 (we are referring, in particular, to the representation in 
the second line of the cited equation), we readily infer 


r 'Ku _ Tr' u i(>) I Tr' u l<) 

ui ,0 — 01,0 ' 01,0 


where 


W-)._ ^ 

Wo •“ a 


/ dx 2 ...dx d d x i y iD^(x,y) 

’(0,a 2 )x...x(0,a d ) 2 


y=x, x 1 =0 


for • e {>, <}. 


(3.91) 


In the next two paragraphs we introduce convergent series expansions for the func¬ 
tions 5i(o > ' ) and ; these expansions ultimately yield the analytic continuations 
of these functions and of at u = 0 (and so, they determine the renormalized 
force on 7r lj0 according to Eq. (3.89)). We also give remainder estimates for these 
series. 

3.8.1 Series expansions and analytic continuations for S r i,o > \ S r i o < ' > 


In¬ 


serting into Eq. (3.91) the series (3.31) for d x i y iD <y s > \yi, y) and integrating term by 
term, we have 


01,0 




2 ai T(^ 


» + i, E 

n6N d ^ 


I Ui 7T 


Y 


«i ) 


w„ (u+1) rf — ,u^T 


(3.92) 


The above expression can be proved to converge for al l u G C (by a simple varia tion 
of the proof of the convergence of the expansion (3.79) for E u, ( > '\ see paragraph C.l 
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of Appendix C); thus, Eq. (3.92) gives the analytic continuation of Sio > ' > to the 
whole complex plane, in particular a t u - 0. 


On the other hand, using the series ( |3.37 ) for d, 
tion (3.91), we can show with some effort that 


x 1 y 1 B s 


(<), 


x, y) along with the defini- 


~«,(<) _ _ 
n i,o — 


. u —3 

k Tw 


E 


■i) 


d—p 


2 i+lr (T) u (rf-p)SCp-(<rr)5 


E 


&cr,p 


E 

hezp L 


((Zl /ll) P u—p—3 


P?,p(h) 


o&Sd 

T ^ 

-Vu_^ 


in the above we have put, for brevity, 

S d '■= {deS'rf | d(l) = 1} , a s> i := 1 , 


P?,p( h ) 


T 


Ba, l(h) := (ai/ii) 2 


(3.93) 


(3.94) 


:= II a *W ’ B °A h ) '■= ( a i^i) 2 + y^( a a(b ^) 2 for creSd, pe {2,..., d} . 


i=2 


i=2 


see 


The last result is derived in a manner similar to expansion (3.80) for E u ’bO ( 
Appendix [Bj; besides, there hold considerations analogous to those below the cited 
equation. More precisely: the terms in the right-hand side of Eq. (3.93) with 
p G {1, ...,d} and h = 0 have a simple pole at u = p + 1; after removing these 
finitely many terms, the series in th e right-hand side of Eq. (3.93) converges for 
all u G C. Summing up, Eq. (3.93) gives the analytic continuation of 5“o < ' ) to a 


meromorphic function on the whole complex plane, with simple poles at 

u G {2, 3,..., d, d + 1} . (3.95) 

Summing up, u = 0 is a regular point for the analytic continuations of both ^ 

and 5io i now, recalling the definition (3.89) for the renormalized total force on 
7r li0 and comparing with Eq. (3.91), we get 


rgren r*u ,(>) 

01,0 — 0 1,0 


«=0 


+ 


u,(<) 


u=0 


(3.96) 


where the two addenda on the right-hand side simply indicate the expressions (3.92) 
and (3.93) evaluated at u — 0. 


3.8.2 Estimates for the series (3.92) (3.93). The arguments of paragraph 
and Appendix 0 for the expansions (|3.T9|) (|3.80[) of E u ^ and E u ^ <s > can also 


3.7.2 


be adapted to deduce remainder estimates on the series (3.92) (3.93) for 5i.o > ’ ) an d 
5io < ' > - Here we only report the final results, concerning truncation at a suitable 


order N ; in our presentation we adopt the notations introduced in subsection [374 
particular, Eq.s (3.45) ( 3.46[ )). 


see, m 
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For N e (0, +oo), let us write 


Tr n i(>) . 

U 1,0,AT ■- 


>«,(>) 


2 ai r(^) ^ 

v 2 ' neN d ,|n|^Af 


7r n i(>) _ 7r u i(>) I o> u i(>) 

u 1,0 — U1,0,JV "r ' /v L,0,JV :) 


1 ’°’ JV ' 2aiT(^) f- 

v 2 ' neN d ,|n|>AT 


E ^ W->r 


M +1 2 rp 

2 ,UJnT 


U + l 2 




(3.97) 


the remainder 93“’^ has the bound 


^ 1 , 0 ,N 


, ( »I /c 8 * u 7t 1- ^ max(a s “ +1 , A Ru+1 ) 


2 rf + 1 a?|r(^±i 


HP <*, 


7T 2 T Uhu+1 
2 


, 1— 


for either ^ 1, N > 2 \fd or 9?w < 1, N > 2Vd, + —\———— 

7r V 2 aT 


(3.98) 


Moreover, let 


7r u >(<) _ '?r u >(<) j_ 

v'1,0 — u i,o, 1V ~r - n l ,0,N ) 


7r n i(<) . 

U 1,0,AT •" 


u —3 

K T^~ 


2 d+1 r(^) ^ (d-p)!(p-l)! ^ (ttT)I 


^ (aihi) 2 P, 


>«.(<) ._ 

v i,o,v - — 


hGZP,|h|^Af L 
k u T^t 


B,A h)\ T 

__ __ P2={=1 


-Bg,p(h) 

T 


''—l i rf -p 


a; J 2 f—ij Uct, p 

2 d+1 r(^±I) ^ (d—p)!(p—1)! (7rT )2 


hGZP, |h|>AT 




B«,(h)\ r 

-2 V ^ 


Bg :P ( h) 

T 


(3.99) 


(recall that |h| := /?, 2 + ... + ); as for the remainder 93“ (containing no 

singular term with h = 0), there holds the estimate 
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K 


5R u 


E 


2 d+i |r(^)| ^ (d—p)\(p—i)\ ^n p/2 


E &a,p 

7rP/2 


a 2 max(a s ^- 3 , A^" 3 ) H™ (a, 


a 2 p+3 — Ku 


, Ku — p— 1 + 


(3.100) 




a 2 p+1 —Km 


T 


,Ku—p—1 


for either Km ^ 2, N > 2\fd or Km > d+ 1, N > 2 y/d H—y---— . 

a V 2a 


Concerning the case 2 < Km ^ d + 1, not taken into account in Eq. (3.100), 
there hold considerations analogous to the ones below Eq. (3.87); also in this case 
the corresponding reminder estimates would involve rather cumbersome expressions 
which we choose not to discuss here in view of the fact that we will be interested 
only in the case u = 0. 


The evaluation of the force on vryo presented in the subsequent subsections 4.1 4.2 


for d — 1, 2 will be based on the truncated expansions and on the related remainder 
bounds discussed in this paragraph. 


3.9 Scaling considerations. From Eq.s (3.55 3.58) and from the expressions for 
the Dirichlet functions given in subsection 3^3[ we easily infer the following relation 
for each component of the stress-energy VEV (/i, v G {0,..., d})\ 


(0|T“(x)|0> = < 


T;(x*;p) , 


(3.101) 


where is a suitable function and x*, p are, respectively, the d-tuple and the 
(d — l)-tuple with compoments 


x 


x* := — G (0,1) for i G {1,..., d} , 


Pi := — for i G {2,..., d} . (3.102) 

CL\ 


For d = 1, the variables pi are not defined and only depends on x\ = 


Similarly, for the regularized pressure acting on any point x in the 
side 7T 10 , we deduce from Eq.s (3.60) (3.101) that 

for i G {1,..., d} 


interior of the 


p“(x) = < d 1 v r( x *;p) 


(3.103) 


11 


It is apparent from Eq.s (3.8) (3.11) that the eigenfunctions Ek(x) := F/ Cl (x 1 )...Fk d (x d ) and 


the corresponding eigenvalues can be written in the form 


n ? ;7T 


F k (x) = cq 2 y n ,p(x*) , w k = a 1 A n p for k, = — (ie {1,..., d}) , 

CLi 

for some suitable functions <p n .p, and some coefficients A n p. Using the eigenfunction expansion 
for the Dirichlet kernel (see Eq. (3.19) in Part I), we obtain 

D s {x,y) = Gq (d ' 2 )s £> Sj p(x*,y*) , 
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where p“ are suitable functions and x* is defined as in Eq. (3.102) at points on the 


boundary. Clearly, the same conclusions can be drawn for the pressure on any other 
side tt P iA (p g {1,..., d}, A G {0,1}). 

Analogous considerations hold for the total energy and integrated force on the 
boundary of the spatial domain. On the one hand, concerning the bulk energy, 


from the expansions derived in subsection 3.7 we easily infer (indicating with E u a 
suitable function) 

E u = a u 1 - 1 E u {p) . (3.104) 


On the other hand, as for the total force on 7r 10 (for example), from Eq.s (3.92) 
(3.93|) it follows that 

ffi,o = ar 2 Fi,o(p) (3-105) 


for some suitable function F“ 0 . Again, similar results hold for the total force on any 
other side 7r P)A . 

By analytic continuation at u — 0, we obtain the renormalized counterparts of the 
above relations: more precisely, we have 


(0|f^(x)|0) ren = a x (d+1) T mv (x*; p) , 
E ren = ah 1 E(p) , 


K en ( x ) = a i (d+1) Pt(x*;p) , 

K o=ar 2 F ll 0 (p) 


(3.106) 


(where the right-hand sides of the above relations are obtained evaluating at u = 0 
the functions in the right-hand sides of Eq.s (3.101 3.105). 


Due to the remarks of this subsection, for any spatial dimension d the analysis of the 
renormalized stress-energy VEV, total energy, pressure and of the integrated force 
can always be reduced to the case oq = 1; we will use this fact in the next section 
on the cases d — 1 and d — 2. 


4 The previous results in spatial dimension d £ {1,2} 


4.1 Case d = 1 (the segment). As a first application of our framework, let us 
consider the case 


d = 1 , ft = (0,ai) (di > 0) . 


(4.1) 


As a matter of fact, we already analysed this configuration in Section 6 of Part I; 
the present Eq. (4.1) corresponds to Eq. (6.1) of Part I (with a = cp). 


for some suitable function 'D s ,p' 1 from the above relation we can easily infer that, for any pair z , w 
of spatial variables, 


dzwDsix, y)=a 1 (d 2)s 2 d ZtWt ® s n(x*, y*) 


Using these results, one can easily deduce Eq. 
subsection. 


(3.101) and the subsequent statements of this 
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In Part I we performed the exact computation for the renormalized VEV of the 
stress-energy tensor and of the pressure for various types of boundary conditions. 
Here we carry out an approximate evaluation of (0|T , /iJ/ (x)|0) re n, Pi en ( x )U!=o an d of 
E ren for Dirichlet boundary conditions, truncating the series expansions for these 
quantities derived in the present work for a box in arbitrary spatial dimension. Our 
aim is just to check the validity of the general methods developed here; to this 
purpose, we compare the results obtained for the renormalized stress-energy VEV, 
pressure and total energy, respectively, with those reported in Eq.s (6.24), (6.26) 
and (6.27) of Part I. 

In our computations we only consider the case with 


a\ = 1 


(4.2) 


which cause no loss of generality due to the scaling considerations discussed in 


subsection 3.9 (of course, due to Eq. (4.2), the rescaled variable x\ := x 1 /'oq in fact 
coincides with x 1 ) ( 12 ). 

Before proceeding to the evaluation of the renormalized VEVs (0|T Mi ,|0) ren and p\ en , 


let us recall the representation (3.56 3.58) of the regularized stress-energy VEV in 
terms of the (>) and (<) parts of the Dirichlet kernel; these parts depend on the 
choice of a parameter T > 0, which however has no effect on the sum D s = D + 
D After fixing T, we can approximate D i > ' > and D truncating their series 
expansions at some sufficiently large order N, giving estimates on the remainders as 
well. Analogous considerations hold for the spatial integral of the diagonal Dirichlet 
kernel, ultimately giving the renormalized bulk energy E ren (see Eq.s (2.36) (3.78) 


and (JA83J)). 
Here we choose 


T = 1 


N = 5 


(4.3) 


the truncated sums are evaluated numerically, and the remainder estimates (J3.49 
3.54) and (3.85) (3.87) are also taken into account, fixing ^ 


a = 0.03 . (4.4) 

Of course, the bounds obtained for the remainders associated to D i > ' 1 and D 
allow us, in turn, to infer error estimates for the approximate expressions of both 
(0|f^|0) reri and pj( en . 

Let us first consider the renormalized stress-energy VEV (0|T Ali/ |0) rer j; according to 


the analysis of subsection 3.5, this is obtained setting u — 0 in Eq.s (3.55 3.58). 


12 Let us stress that, for ai = 1, the quantities (0|7k„|0)r 

and pl en (see Eq. (0061)). 


the rescaled functions TJ)®" 

13 We make the choice (4.4|) for a because it is close to the value minimizing the error esti¬ 


and p\ en are, respectively, equal to 


mates cited above; a more precise evaluation of this optimal value for a would require a laborious 
numerical analysis which we prefer to avoid here. 
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In reporting our results, we distinguish between the conformal and non-conformal 
parts of each component, which are respectively denoted as usua l wit 
and (0\T&\0) 


i (0|T^ } |0) 


(see Eq. (2.12)); notice that (since d — 1) Eq. ( 2.13[ ) gives 

6-o. 


(4.5) 


The graphs of the conformal and noncomformal parts of each stress-energy com¬ 
ponent are shown in Fig.s [I] and [2] (recall that we refer to the rescaled variable 
x\ := x 1 ja\ = a; 1 ). Let us comment briefly on the above graphs. Apart from 



Figure 1: d= 1: graphs of (0|f 0 ( Q } |0) ren and (0|T^" ) 10)^^. ■ 
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Figure 2: d= 1: graphs of (0|T 1 l 6|0)ren and (OlT^IO)^ ■ 

(0|Tg"^|0) ren , it appears that all the other components of the renormalized stress- 
energy VEV are constants and (0|T^ ! |0) ren is very small; indeed, our computations 
with N = 5 (T = 1, a = 0.03) ensure, for all x = x\ G (0,1), 


(0|f o ( o O) (x)|0) ren = (0|T^(x)|0) ren = -0.1308997±8- 10 


-( 0 ) 


-7 


|(0|f 1 ?)(x)|0) ren |^4.10 


-6 


(4.6) 
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Concerning (0|To"^(x)|0) ren we have, for example, 
(0|T 0 ( 0 ' > (x)|0) ren = 


3.141593 ± 10~ 6 for = 1/4 
1.570796 ± 10~ 6 for x\ = 1/2 


(4.7) 


These results are in agreement with the exact calculations of subsection 6.6 of Part 
I, which gave the following outcomes (see Eq. (6.24) of the cited subsection with 
a = cq = 1) ( 14 ) 


{0|f 0 ( „ O, (x)|0)„„ = (0|T™(x)|0)„ n = ~ , 




7r 


(0|T o ( o ' > (x)|0) ren — 


7T 


2sin 2 (7rx 1 ) 
Next, let us pass to the evaluation of 

^en(O) :=PrW 


(0|f 1 ? ) (x)|0U = 0. 


(4.8) 


34=0 


(4.9) 


this is nominally the “pressure” on the boundary point x 1 = 0 but in fact coincides 
with the force on this point, due to the zero dimensionality of the boundary. 

Let us consider the prescriptions (3.62); computing the derivative of the functions 

obtain 

(4.10) 


D[% Dyl appearing therein with the choices a\ — 1, T — 1 and N = 5, we obtain 


F ren (0) = 0.1308997 ± 3 • 10 


-7 


(again, the error is obtained using the remainder estimates of subsection 3.8.2 with 
a = 0.03). The above result is in agreement with the exact expression F ren ( 0) = ^ 
derived in subsection 6.6 of Part I (see Eq.s (6.24) (6.27), and set a = a\ — 1 
therein). 


Finally, we consider the renormalized bulk energy E ren ; the series expansions (3.79) 
(3.80) derived in subsection 3.7 (with the previous choices of ai,T,N) allow us to 
infer 

(4.11) 


-13 


E ren = -0.1308996938996 ±3-10 

(where the error is obtained using the remainder estimates of subsection 3.7.2[ again 
with a = 0.03). The result (4.11) agrees with the exact computation E ren = —^ 
obtained in subsection 6.6 of Part 1 (see Eq. (6.26), and set a = a\ — 1 therein). 


14 


To make a comparison with Eq.s (4.6) (4.7), note that 


-= -0.13089969... , 

24 


2shr(7ra; 1 ) 


3.14159265.. . for i 1 = 1/4 

1.57079633.. . for x 1 = 1/2 
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Let us stress that our approximants by truncation, converge quite rapidly to the 
exact results; in order to exemplify this statement, we notice that, by slightly in¬ 
creasing the value of the truncation order N, we obtain a remarkable improvement 


of the error estimates. For example (using again the estimates of subsection 3.7.2 
with T — 1 and a = 0.03) the error ±3 • 10~ 13 in Eq. (4.11) (for N = 5) becomes 


±2 • 10" 46 for N = 10, 6 • 1(T 101 for IV = 15 and 4 ■ lO" 177 for N = 20. 
4.2 Case d — 2. Let us now pass to the 2-dimensional case: 

d — 2 , n = (0,ai) x (0, a 2 ) (aq, a 2 > 0) . 


As in the previous subsection, we fix 


a i — 1 


(4.12) 


(4.13) 


and consider different values of a 2 ; let us repeat that the above choice does not 


imply a loss of generality, due to the scaling properties of subsection |3.9| Moreover, 
we present the final results in terms of the rescaled coordinates x\ := x l /a\ = x 1 , 
xl := x 2 /a 2 E (0,1), defined in Eq. (3.102) ^ 


Again, the basic elements to compute the renormalized stress-energy VEV and the 
pressure are the Dirichlet functions D^, along with their spatial derivatives, 


for which we use the truncated expansions (3.49 3.54) and the remainder bounds of 


Eq.s (3.49 3.54). 


Needless to say, analogous considerations also hold for the renormalized bulk energy 


and for the integrated boundary forces (see subsections 3.7 and 3.8, respectively). 


Let us first consider the stress-energy VEV and the pressure; as examples, we com¬ 
pute these observables for the two configurations with 


a 2 — 1 


and 


a 2 = 5 . 


(4.14) 


In these cases, for the parameter T of the decomposition into (>) and (<) parts and 
for the truncation order N, we make the following choices: 


T = 1 
T = 1 


N = 7 
N = 9 


for a 2 — 1 ; 
for a 2 = 5 . 


(4.15) 


therein the choice ( 16 ) 


The truncation errors in Eq.s (|3.49j)3.54|) are evaluated making for the parameter a 

(4.16) 


a = 0.04 


15 

n ren 

Pi 


Similarly to what we said in the footnote 


12 


on page 


35 


for ai = 1, the quantities (0|T M „|0) 


rem 

Ten 


and 37™ (to be discussed hereafter) do in fact coincide with the rescaled analogues T, 
pf en , FJ™ introduced in Eq.s (3.101) (3.103) (3.105)). Besides, the lenght a 2 of the second side is 
identified with the ratio p 2 (see Eq. (3.102[)). 


16 This choice can be justified by considerations similar to the ones in the footnote 13 of page 35 
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The renormalized stress-energy VEV (0|T A1 ^|0) reri is obtained setting u = 0 in Eq.s 
(3.55 3.58). Again, we separate the conformal and nonconformal parts, respectively 


indicated by the superscripts (0) and 
dimensional case, 

6 = 1/8 . 


; recall that Eq. (2.13) gives, in the two- 


(4.17) 


In the following we present the graphs for (0|T / ( //|0) ren and (0|7~^™ ; |0 ) ren obtained 
from the previous truncated expansions; more precisely, Fig.s[3}j5]and Fig. s[6](9] show, 
respectively, the results obtained for the configurations with a 2 = 1 and a 2 = 5. 
In the cited figures we refer to the variables x\ := x l /a* € (0,1) and, keeping into 


;(■)! 


account some obvious symmetry considerations ( i( ), we only show the graphs for 


4 G (0,1/2) fori G{1,2}; 


(4.18) 


moreover, in the case of a square box with oq = a 2 = 1 we do not report the graphs 
for the conformal and non-conformal parts of (0|T 22 (x)|0) reTl , since these are equal 
to the corresponding parts of (0|Tn(x)|0) ren . 



Figure 3: d = 2: graphs of (0|f 0 ( o } (x)|0) reri and (0|T o ( " ) (x)|0) re „ for a 2 = 1. 


1 ' Indeed, every component of the stress-energy VEV can be shown to be symmetric under the 
exchange x l a* — x l (or x\ 1 — x\) for i G {1, 2} . 
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Figure 4: d = 2: graphs of (OlT^f^x)^)^ and (0|T^"^(x)|0) rera for a 2 = 1. 




Figure 6: d = 2: graphs of (0|TQQ^(x)|0) ren and (0|To"^(x)|0) rer i for a 2 = 5. 
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<0|jft|0 > 



Figure 7: d = 2: graphs of (OlT^^x^O)^ and (0|T^"^(x)|0) r en f° r «2 = 5. 





Figure 9: d = 2: graphs of (OjT^^x^O^en and (0|T2"^(x)|0) re n for 02 = 5. 
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Concerning the error estimates, for /i, v G {0,1, 2} and • G {0, ■}, let us introduce 
the following notation: 



remainder corresponding to our approximation 
by truncation of (01 10) T . eT1 . 


(4.19) 


For a 2 = 1, our choice N — 7 (T — 1, a — 0.04) yields the uniform bounds 


|4 0) I<2-io 


-12 


If, 


| 4 0) | < 6-10 


-12 


1^ 


(■)| 


V 


<3-10 


>K2- 10” 11 ; 

for i,j G {1,2} 


For a 2 = 5, our choice N — 9 (T — 1, a — 0.04) ensures 


( 0 ), 

oo I 


If, 

144 ^2-io 


<5-10 


-12 


If, 


(■)| 
00 I 


<4-10 


-li 


-li 


If 


(■)i 
ij I 


< 8-10 


-li 


for i,j G {1,2} 


(4.20) 


(4.21) 


Now, let us evaluate the pressure p[ en (x) at points x of one side; as in the construc¬ 
tion of the general theory we consider, as an example, the points x = (0,x 2 ) in the 
interior of the side 711 , 0 , making reference to the prescription (3.62). 

Fig. 10 shows the graphs obtained for p} en (x) as a function of x 2 := x 2 /a 2 (again, 
choosing T — 1 and truncating the related expansions to order N — 7, for a 2 = 1, 
and N = 9, for a 2 = 5). 



Figure 10: d = 2: graphs of p[ en for a 2 = 1 (left) and a 2 = 5 (right). 


As for the error, indicating with e± the remainder associated to our approximation 
by truncation of p [ en , we obtain the following uniform estimates (setting a = 0.04): 


ei | ^ 2 • 10 12 for a 2 = 1 ; 

ei|^5-10' 12 for a 2 = 5 . 


(4.22) 


42 














Before moving on, let us briefly comment on the behaviour of the renormalized 
pressure p\ en near the edge x = 0. Indeed, specializing to the present d = 2 case 


the considerations of paragraph 3.6.2 one can prove that, for all a 2 > 0, 
1 


K en (x) = 


+ 0((x 2 ) 2 ) for x = (0, x 2 ) and x 2 —> 0 + . 


(4.23) 


327r(x 2 ) 3 

Now, let us pass to the computation of the bulk energy and of the integrated force. 
For each one of these two observables we consider several configurations, corre¬ 
sponding to different values of a 2 ; for any one of these values, we consider the 
decomposition into (>) and (<) parts, and choose the truncation order N of the 
related expansions so as to obtain error estimates all of approximatively the same 
order of magnitude, fixing again 

a = 0.04 . (4.24) 

In order to obtain the renormalized bulk energy E ren , we first consider its reg¬ 


ularized version (3.78), along with the series expansions (3.79) (3.80), giving the 


analytic continuations of the functions E u ’*•>) and E u 4<), respectively. Due to the 


considerations of subsection |3.7[ we can simply put u — 0 in these expansions, since 
no singularity arises there; next, we truncate the corresponding series at a suitable 


order N (see the comments above) and use the remainder estimates (3.85) (3.87). 
In this way we obtain, for example, the following results: for ri\ — 1 (and with the 
choice T = 1, as usual) 


E ren = _ 1 .73 691 776± 4 - 10 -8 for a 2 = 0.1 (N = 40) 
E ren = + q.03524178 ± 3 • 1(T 8 for a 2 = 0.5 (N = 8) ; 

E ren = + q.04104060 ± 9 ■ 1(T 8 for a 2 = 1 (N = 4) ; 

E ren = _ q.05412096 ± 2 ■ 10~ 8 for a 2 = 5 (N = 7) ; 

E ren = _ 0 .17369178 ± 1 ■ 1(T 8 for a 2 = 10 (N = 15) 


(4.25) 


Again for a\ — 1, one can plot the renormalized bulk energy E ren as a function 
of a 2 ; the graph in Fig. [Tl] has been obtained using the truncation order N = 50 
(with T — 1; with the choices made for the parameters N,T,a\, we know that the 
remainder is smaller than 2 • 10 -3 for a 2 G [0.05,10)). 

Let us discuss some facts regarding the function a 2 i—> E ren (a 2 ), which can be read 


from the graph in Fig. 11 (the results reported hereafter are obtained using standard 
numerical methods, implemented in Matheraatica) ( 18 ). 


18 The derivation of the results in items i)-iv), including the error estimates, should be accounted 
for; as an example, let us give some details on Eq. (4.26). To obtain this result we have used for 


E ren an approximant by truncation at order TV = 50 that, according to the reminder estimates 
(3.85) and ( |3.87 l, gives E ren up to an error ^ 10~ 272 for a\ = 1 and 0.5 < a 2 < 1. This large order 
approximant of E ren has been maximized numerically with respect to a 2 via Mathematica, asking 
for a precision of order 10“ 15 on the maximum point. The final results have been prudentially 


truncated to 8 digits, yielding Eq. (4.261; in view of the previous considerations, they are very 
likely even though not certified. 
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Figure 11: d = 2: graph of E ren as function of <22 . 


i) There is only one point of maximum a™ ax ; our approximation by truncation at 
order IV = 50 gives 


a max = 0.72719110 ± 10“ 8 , 

E ren (a™ ax ) = 0.04472675 ± 10" 8 . 


(4.26) 


ii) E ren vanishes for two values op 1 < ap of 02 ; these are found to be 

c4 1} = 0.36538151 ± 10~ 8 , 
op = 2.73686534 ± 10" 8 . 


(4.27) 


E ren is positive for aP < a 2 < aP and negative elsewhere. This feature was also 
pointed out in [22]; therein it is stated that op = (a^) -1 , a relation (approximately) 


verified by the numerical values in Eq. (4.27). 

iii) For <22 —> 0 + , E ren has the asymptotic behaviour 

E ren (a 2 ) = 7 ~~~“T 7 r (l + 0(a 2 )) with e 0 = -0.02391 ± 10 

(a 2 ) 2 V ) 


(4.28) 


iv) There are indications that E ren approaches an asymptote for a 2 —» + 00 ; taking 
into account values of the abscissa up to a 2 = 100, we find that this asymptote is 
the straight line 


y = m,E a 2 + qE with 


f m E = - 0.02391416 ± 10" 8 
\q E = + 0.06544985 ± 10" 8 


(4.29) 


Finally, let discuss the renormalized total force 3™ acting on the side 7Ti i0 ; following 
the analysis of subsection 3.8, we consider the expression (3.96) and represent the 


functions appearing therein using the series expansions (3.92) |T93). 

$[ e Q is obtained setting u — 0 in these series expansions; for the actual calculation, 
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these can be truncated at a fixed, sufficiently large order N using the remainder 


estimates (3.98) (3.100). 


In this way we obtain, for example, the following results: for aq — 1 (and with the 
choice T = 1, as usual) 


jren = + 2.3914162 ± 4 • 10“ 7 

= + 0.0956400 ± 2 • 10~ 7 
= + 0.020520 ± 6 • 10 -6 
= - 0.173692 ± 4 • 10" 6 
= - 0.412833 ± 3 • 10 -6 


for a 2 — 0.1 
for a 2 = 0.5 
for a 2 = 1 
for a 2 = 5 
for a 2 = 10 


(N = 50) ; 
(N = 9) ; 
(N = 4) ; 
(N = 7) ; 
(N = 15) . 


(4.30) 


Again for aq = 1, one can plot the renormalized force as a function of a 2 ; the 
graph in Fig. 12 has been obtained using the truncation order N = 60 (with T — 1; 


with the choices made for the parameters N, T, op, we know that the remainder is 
smaller than 2 • 10~ 3 for a 2 G [0.07,10)). 


F'~ 



Figure 12: d — 2: graph of as function of a 2 . 


In the following we briefly discuss a number of facts concerning the function a 2 t —> 
$[ e o( a 2 ), which appear in the above graph. 

i) The function under analysis is strictly decreasing for a 2 € (0, +oo). 

ii) There is only one value of a 2 , which we indicate with a 2 , where vanishes; 
our approximation by truncation at order N = 60 allows us to infer 

a 2 = 1.3751543 ± 10" 7 . (4.31) 


iii) For a 2 —» 0 + , has the asymptotic behaviour 


iTT 0 > 2 ) = (l + 0(a 2 )) with / 0 = 0.02391 ± 10 


-5 


(4.32) 


iv) It appears that approaches an asymptote for a 2 —» +oo; considering again 
values of the abscissa up to a 2 = 100, the equation of this asymptote is found to be 


y = ma 2 + q 


f m = - 0.0478283 ± 10" 7 
\q = + 0.0654498 ± 10" 7 


(4.33) 
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Let us mention that, in agreement with the general considerations of subsection 4.4 
in Part I, the above results for the renormalized total force on the side 7Ti o of the 
box could be equivalently derived by differentiating the renormalized total energy 
E ren with respect to the lcnght 02 of the edge of the box perpendicular to 7r 10 . 


To conclude, let us compare the previous results about E ren ,$ To* with the calcula¬ 
tions of Bordag et al. [5]. The authors of [5] derive (both by Abel-Plana formula 
and by zeta regularization) series expansions different from ours for the renormalized 
bulk energy and for the force on one side; besides, they give no remainder estimates 
for these expansions ( 19 ). 

The numerical values of E ren , ^\ e Q given by our previous analysis are in good agree¬ 
ment with those arising from the expansions in [S], a fact strongly indicating the 
equivalence between our approach and [5]. Let us also mention that the results of 
[5] about E ren or are equivalent to the ones of |3| HU flfi j. 


Acknowledgments. This work was partly supported by INdAM, INFN and by 
MIUR, PRIN 2010 Research Project “Geometric and analytic theory of Hamiltonian 
systems in finite and infinite dimensions”. 


19 In fact, a not so trivial analysis (which we do not report here for brevity) allows to conclude 
that the series expansions given in [5] also converge with exponential speed. 
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A Appendix. Absolute convergence and remain¬ 
der bounds for the series of subsections 


3.3 


Let us consider the general f rame work of the subsections mentioned in the title. 
Note that the representation (3.30) for the function Z^i > ' ) (x, y) contains the series 


neN d 


Sr 0> W n T ) Cn(*,y) 


(A.l) 


On the other hand, in the expansion (3.36) for D^^x, y) there appears a series 


over hGZ“ and 1 g{ 1,2} , that we have reexpressed in the equivalent form (3.39) 


We pointed out in paragraph 3.3.2 that the terms in this series with &hi( x , y) = 0 


become singular at s = d/2; these have either |h| = 0 or |h| = 1 (see Eq. (3.22)) 
and, after removing all terms with |h| = 0,1 we are left with the series 


E s ' v - 


hez d , |h|> 1 
ie{i,2} d 


M x ,y) 

T 


(A.2) 


Note that the absolute convergence of the series (A.2) is equivalent to the absolute 


convergence of the series in the second line of Eq. (3.39) 


In the next subsections we prove the convergence of the series (A.l) (A.2) for all 
s 6 C , deriving as well remainder estimates for b oth o f them ; the statements in Eq.s 


( |3.30| ( |3.50| (related to D^ipc, y)) and in Eq.s ( |3.36| ) ( 3.53[ ) (related to D^^x, y)) 
follow easily from these results. 

The arguments presented in this appendix can be easily generalized to derive the 
analogous conclusions for the series in Eq.s (3.31) ( |3.51 ) and (3.37) (3.54) related to 
the derivatives d zw D^ {x., y) and d zw Di < \x., y); we will briefly return on this topic 


in the sequel (see subsections A.2 and A.3) 


In the final subsection of this appendix we also derive the asymptotic (3.47) for the 


function Hffl in our remainder estimates. 


Before proceeding, let us recall the following definitions (see Eq. (3.45)): 


a : = min {m|, A := max jail, Iz1 2 := > zf for z = nGN^ or z = hGZA 

1=1 

A.l Preliminary estimates. First of all, let us consider the upper incomplete 
gamma function (see Eq. (3.29)) 

r+oo 

F(s,z):— / dw e~ w w s ~ x (s G C, z G (0, +oo)) 
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and note that the above definition implies the following: 


|r(s,z 2 )| < r(9fa,Z!) for all sE C, zi, z 2 E (0, +oo) with 0< z\ ; (A.3) 

|T(s, z) | ^ (1 —a)^ s e~ az r(5Rs, (1 —a)z) for all sE C, zE (0, +oo), aE (0,1). (A.4) 
Hereafter we are going to prove that 


E |h| p |r(s,/3|h| 2 )| ^H$(a,P;$ls,p) 

heZ d , |h|>AT 

for all s E C, N, p E R, /3 E (0, +oo), a E (0,1) such that (A-5) 

either p ^ 0 and N > 2\fd or p > 0 and N > 2 Vd + J- 

V' 


p 


2a(3 ’ 


where Hffi is defined according to Eq. (3.46), i.e. (for cr E R), 

/3; a, p) : = 


7r 


d/2 


N-Vd 


1 


r(a, {l-a)pN 2 ) V 


f^;a/?(iV-2v/d)d 


V 2 


(l-a) CT (a/3)^r(|) \N- 2 Vd / 

In the sequel we always assume s E C, p E R, NE (2 y/d, +oo), /3e (0, +oo), aE (0,1). 
In order to derive the estimate (A.5) we first point out that, due to Eq.s (A.3) (JA.4 ), 


there holds 


|r(s,/3|h| 2 )| < (l-a)- Rs e~“^ |h|2 r(K5,(l-a)/3IV 2 ) 

for all hGZ d with |h| > AA ; 


(A.6) 


this in turn implies 

E |h|Hr(s,/3|h| 2 )|^ (l-a)" !Rs r(3«s,(l-a) y 0IV 2 ) E \h\ p e~ am2 . (A.7) 

heZ d , |h|>JV h£Z d , |h|>JV 


To proceed, let us note that, for any non-increasing function F : [N — 2\/d, +oo) —> 
[0, H-oo), there holds the following estimate (see, e.g., [28], P- 691): 


E F <i h i)^ 

heZ d , |h|>JV 


2vr d / 2 

W 



(A.8) 


noting that r+Vd = + ^ or r > N — 2V~d > 0), the above 

inequality yields 


vd—1 


E F <i h i)^ 


r»+oo 


hez d , |h|>JV 


2t r d / 2 ( N-Vd 
r(|) liv— 2 ^/ Jn~ 2Vd 


dr T d ~ 1 F{r) 


(A.9) 


48 















We want to employ the above bound to estimate the sum in the right-hand side of 
Eq. (A.7); this case involves the function 


F(t) = t p e~ afir \ (A.10) 

which is decreasing on [N — 2yfd, +oo) for 

either p ^ 0 and N> 2 y/d (as before) or p > 0 and N> 2 \[d + , / —. (A. 11) 

V 2ctp 

In each one of the above two situations, making the chang e of variable v := a f3 r 2 G 
[a f3(N — 2\fd) 2 , +oo) and recalling the definition (3.29) of the upper incomplete 
gamma function, we easily obtain 


f+oa -i 

/ dr r d+p ~ l e~ a ^ r2 = -^ 

’ N—2Vd 2 (a (3) 2 


^; Q| 9(lV-2v / rf) 2 ') , (A.12) 


which, along with Eq.s (A.7) (A.9) and (A.10), proves Eq. (A.5). 


In passing we point out that Eq. (A.5) implies the analogous estimate (P 


l n ! P l r ( S ’^l n ! 2 )l ^ Y d H N\ a i p) 

ngN d , |n|>JV 

for all s G C, N, p G R, (3 G (0, +oo), a G (0,1) such that 
either p ^ 0 and N > 2 \fd or p > 0 and N > 2 \[d + ^ 


(A.13) 


P 


2ap 


A.2 The series (A.l); connections with the expansions (3.30) (3.31) for 
and its derivatives. For sGC and N G (0, +oo), let us define 


77 


(>) i 

s.N 


x,y):= UJ n 2sr ( s ^i T ) C n (x.,y) 

nGN d , |n|>JV 


(A.14) 


of course, the series (A.l) is absolutely convergent if and only if 77^(x, y) < +cx) 
for some N . The rest of this subsection will be mostly dedicated to evaluating 


^S( x >y) • 


■'s.AT 


First of all, note that from the definitions (3.11) of C n ,u4 and (3.45) of a, A it easily 
follows, for all n G N cZ , 

|Cn(x,y)| ^ 1 , (A.15) 


20 Indeed, it sufficies to observe that for any function F : (0, +oo) —► [0, +oo) there holds 

E E F (H) • 

hez d ,|h|>iv nSN d ,|n|>iV 
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(A.16) 


7T 2 2 2 7T 2 |2 

^ M> < <4 < ? M* 


Using Eq. (A.16) we infer, for all s G C, n G N rf , 

2 max(a 2#s , A 2SRs ) 

— 2 — ! -- 


7T' 


■25fts 


n 


-25Rs 


|r( s ,^r)| (Ks, 


7T 2 T 

“A 2 " 


n 


(A.17) 


(to deduce the second inequality we also employ Eq. (A.3)). Returning to the 
definition (A. 14), we obtain 


"SfcrX E iW*. 


neN d , |n|>JV 


7T 2 T 

~A^~ 


n 


(A.18) 


and using the estimate (A.13) we conclude the following, for any a G (0,1): 


( » M / TT 

/c s,jva x > y) ^ 2 d 7r 25Rs Un ' ’ - 


T 


a, ; 3?s, —23fts 


(A.19) 


for either U?s ^ 0, A - > 2 \fd or < 0, N > 2 y/d + — \ . 

7T V aT 

As noted before, the finiteness of 77^(x, y) proves the absolute convergence for all 
s G C of the series (A. 1), contained in E q. (3.30) for Di > \x, y); the bound ( 3.50 ) 
for the remainder R^(x, y) of Eq. (3.49) follows straightforwardly from Eq. (A.19) 
since we have 

|fl S (x ' y)l " ai ...a!|r( s )i R $ (x ’ y) ■ ( A - 2 °) 

Let us mention that the series in Eq. (3.31) for d zw D^ can be discussed similarly; 
in place of 77^(x, y), we use the quantity 

(*,u>)ft£jv(x,y) := w n 2s r (s,ulT) d zw C n (x,y) 

neN d ,|n|>W 

and we estimate it similarly to what we did with 77^(x, y), noting that 

7T 2 


(A.21) 


|cW? n (x,y)| ^ |n| 2 


(A.22) 


The conclusions of this analysis are the absolute convergence for all s G C of the 
series (3.31) and the remainder bound (3.51). 
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A.3 The series ( A.2[): connections with the expansions (3.36) (3.37) for 


and its derivatives. Keeping in mind the definitions (3.16), (3.18 


(3.33) of b hh 5i and V s (and noting that |^| = 1), let us put 

M x > y) 


and 


^S( x >y) : = X] 

hez d , |h|>v 
ie{i,2} d 


V s _i 
S 2 


T 


(A.23) 


for sGC and IV G (0, +oo). Let us consider the series (A.2): clearly, this converges 
absolutely if and only if 7?.^(x, y) < +oo for some N. 

Most of the sequel will be dedicated to evaluating 7Z^(x,y). To this purpose we 


first note that, for 1 G {1, 2} rf , h G Z d and |h| > N > y/d, ( 21 ) 

a 2 f 1 - l h | 2 < &hi(x,y) ^ A 2 fl + ^ 

Using these bounds we easily infer, for all s G C, 

|6hi(x,y) s KC'^(K S ,A) \h\ ms 


(A.24) 


(A.25) 


where, as in Eq. (3.48), we have put (for er G R 
Ci d) A {cr, N ) := max 


•i-flUT#'" 


Next, notice that the identity (3.34) implies 

\Vs(P)\ < |/?| s **r(— 3fts,/3) for all s G C, /3 G (0, +oo) ; 


(A.26) 


21 For example, in order to prove the first inequality in Eq. ( |A.24 ), recall the definition (3.18) of 
&hi; moreover, let us put C/fix, y) := (U^ (x l , 2/*))i=i,...,d and note that 

d d 

|Ui(x,y)| 2 := (U u (x\y 1 )) 2 < Y 1 = d • 

2=1 2=1 

Then, for all i G {1,d}, we have 

d 

&ki(x,y) = Y a *( h * ~ u ui x \y 1 )) 2 > a 2 |h- t/i(x,y)| 2 ^ 

2=1 

! (l h l - |Ui(x,y)|) ^a 2 (jh|-\/d) ^ a 2 ^1 - |h| 


^ a z 
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in particular, from the above relation and Eq.s (A.3) (A.24) (A.25), we deduce 


VA b f 


< 


|25Rs 




r -g?s, 


a 2 (l-f) 


(A.27) 


Inserting the previous results into the dehnition (A.23) of 77^(x, y) we get ( 22 ) 
2“C«(Sft S -f,A0 


f x vl < _ 2,±y J V Ih 

/c s,jvi x >yi ^ jiRs-d/2 2-u i 11 

heZ d ,|h|>iV 

now, using Eq. ( A.13|) we conclude the following, for any s G C and any a G (0,1): 


25Rs -d r ( d.s. /a |h|2 |. (A28) 


7? (<) f x vl < EEETEZ_ 

/c s ,Arl x >yJ ^ T Ks-d/2 


-f, JV) 


( a, 


T 


Vd \2 

N ’ ;|-Ks,2»s-(i 


« 2 (l-f) 2 d 


d d 1 l($ls — -)T 

for either JDs ^ -, A - > 2\/d or > -, TV > 3^ H— 4 ' 2 

2 2' a 


(A.29) 


a 


As anticipated above, the finiteness of 7Z^(x, y) implies the absolute convergence 
of the series ( A.2 ), which appears in the expansion ( 3.36 ) for D^(x, y); besides, 
the bound (3.53) for the remainder R^(x, y) in Eq. (3.52) follows easily from Eq. 
(j A. 29[) noting that 


AA.y)l < (47r )rf/2|r( a ) K*( x -y> 


T#t *-1 


>(<), 


(A.30) 


A similar analysis can be developed for the series in the right-hand side of Eq. (3.37) 
for d zw D^] namely, in place of 7Z^(x,y) we consider 

(,,,»)RS( x .y) := 


E 


heZ d , |h|>AT 
ie{i,2} d 


v -i-Af 


hi \ d z b hi d w b h i 


J*2 


-Vs-i-l ^ 


hi \ '-’zw u hi 


&Z 1 llbb 


T 7 T 


(x,y) 


(A.31) 


and estimate this quantity using Eq. (A.27), along with the following relations: 

\d z ib hl (x,y)\ < A ^1 + |h| , \d ZJU Jj hl {x,y)\ = ^ S. %j ^ ^ 

for z, w = x or y, i,j G {1,..., d} . 

The final result proves, in this case, the absolute convergence of the series in the 
right-hand side of Eq. (3.37) as well as the remainder estimate (3.54). 


22 Note that, since the estimate obtained no longer depends on 1, the sum over 1 £ {l,2} d just 
yields a multiplicative factor 2 d . 
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A.4 Asymptotics for H^\ The upper incomplete gamma function is known to 
possess the following asymptotic expansion (see [29], p.179, Eq.8.11.2): 

T(s, 2 ) = e~ z z s ~ 1 (l + 0(z~ 1 )) for z E R, z —* + 00 , and all sgC. (A.33) 


Comparing with the definition (3.46) of Hffi (for any a E R), we readily infer 

e ~/3N(N-AapVd)jy2cr+p+d-A ^ _|_ 0(A _1 )) 




-fld/2 0<r—2g—Aa (Id 


a(l — a) T(|) 


for N —> +cx) , 


which is the result stated in Eq. (3.47). 

B Appendix. The expansions of subsection 
for the bulk energy 


3.7 


Consider the general framework of subsection |3.7[ where we discuss the regularized 
total energy that, due to the vanishing of the boundary contributions, coincides with 
the regularized bulk energy E u . For the latter we have (see Eq. (3.78)) 

E u = £“’ (>) + E< <] where 


£“>(•).— *1 L 
2 


/ dx 1 ... dx d R^i(x, x) for •€{>,<}. 

'(0,ai)x...x(0,a d ) 2 


Herea fter we show in several steps how to obtain the series expansions (3.79) and 
(3.80) for E u and £7“’^, respectively ( 23 ). 

B.l Derivation of the expansion (3.79). Let us first consider the function 

E u, (>); using the expression (3.30) for the diagonal Dirichlet kernel ■Pj_ 1 (x, x) ap- 

1 2 
pearing in the definition (3.78), we readily obtain 

£«>(>) — 


2 d ~i K n 


ai ...a d r(V) 


IT E <" r 


neN d 


u - 


> W n T 


/ dx 1 ... dx d C n (x, x) . 

' (0,ai)x...x(0,a d ) 


(B.l) 


Recalling the definition (3.11) of C n (x, y), the integral in the above equation can be 
straightforwardly evaluated to yield 


(0,ai)x...x (0,a d ) 


dx 1 ...dx d C n (x, x) = 


i=i Jo 


dx 1 sin 2 ( x l ) = 

c i 


a\ ...ad 
2 d 


(B.2) 


which, along with Eq. (B.l), allows us to infer Eq. (3.79). 


23 Concerning the interchange of certain sums with integrals or derivatives, recall the comments 
of footnote [7] on page 


15 
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B.2 Derivation of the expansion (3.80). Let us now discuss the term 

to this purpose, we consider the expression (3.36) for the Dirichlet kernel D,Z\ (x, y). 

I f 

When evaluated along the diagonal y = x, this expression reduces to 

'bu i(x 


L (x,y) 


rri u - d -1 

i 2 


y=x 


(4vr) d / 2 r(^) 


y, ^ V u-d-i 


hez d ,ie{i, 2 } d 


(B.3) 


here and in the remainder of this appendix, for all h G Z d , 1 G {1,2} d , we use the 
short-hand notation (compare with Eq. (3.18)) 

d 

b h i(x) := 6 h i(x,x) = 'Y^a 2 i {h i - U h {x\x 1 )) 2 . (B.4) 


i— 1 


Eq.s (3.78) (B.3) allow us to infer for E U J < J the representation 

/?“.(<) _ K _ _ n( d ’ T ) 

2(47r) rf / 2 r(^7) ^ 


(B.5) 


where, for suitable sG C and for all T G (0, +oo), we have introduced the function 


Q 


( d,T) 


E * 

hez d ,ie{i, 2 } d ' 


j dx 1 ...dx d V s 

(0,ai)x...x(0,a d ) 


fchi(x) 

T 


(B.6) 


In the next subsection we show that this function can be expressed as 

'B^ p ( h) 


Q (d,T) = y (~V^) d P a Y V 
^ s Z^ ( d-p)\p\ ' *+ 


p =o 


d—p 


cr&Sd hezp 


T 


(B.7) 


where the coefficients a CTi p, B a:P (h) are as in Eq. (3.81). Once (B.7) is proved, this 
equation and (B.5) give the thesis (3.80). 


B.3 Concluding the previous argument: derivation of Eq. (B.7). First of 


all, notice that the definitions (B.6) and (3.33) of Q (d " Tl and V s give 

=[ dr r 3 - 1 Y 5i [ 

JO . i^r. 


dx 1 ... dx d e r h r ) ; (B.8) 


hez d ,ie{i, 2 } d 


'(0,ai)x...x(0,a d ) 


this implies, recalling the definition (B.4) of &hi( x ), 


Q‘* T) = / dr r - 1 J] E S ‘> 


dx i e - ^ {hi-Uliix 1 ,x 1 )) 2 


(B.9) 
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Let us now focus on the general term in the product over i G {1, d}; explicitating 
the sum over k G {1,2} and noting that Ui(x\x l ) = 0, U 2 (x l ,x l ) = x l /di (see Eq. 
d3l6|), we get 


n CLi 


k dx l e -7x( h i- Ul i( xi ’ xi )) 2 = 

hi£Z,li£{l,2} 

( a i h i ) 


E 

hi£ Z 


dx‘e —- / dx‘e~~ 


(B.10) 


The Erst integral in the square brackets above is trivial, since the integrand function 
is constant; moreover 


/ dx l e 


-yr (aihi-x 1 ) 2 _ 


hieZ 


E 


^ CLi hi 


^+00 


(B.ll) 


dz l e ~tt = l dz 1 e ~tt = \ArrT , 


hi£Z J-OO 

where in the second passage we have performed the change of variable z 1 := ctihi—x 1 
Summing up, Eq. (B.10) yields 


f*(Li 


Y 5i > / dx 

hi£Z,li£{l,2} 


i e ~^T (hi-Ui^x 1 ^)) 2 _ a 


e~kk~ ) _ VtttT . (B.12) 


. hiE.Z 


In the following we use the notations introduced in Eq. (3.81). Let us return to Eq. 
(B.9) and consider the product over i G {l,...,d} therein; Eq. (B.12) allow us to 
infer ( 24 ) 


n e *./'«* 


e -^r{h i -Ui i {x\x ' 1 )) 2 _ 


i =1 hi£ Z, li£{ 1,2} 


n 

i =1 


( a i h i) 


di\ 


E 


{-yfkrT) d -P 


Y a -x Y 


Y e tt 

, z 

Z> _ Pf 


— VkVf 




(B.13) 


24 This result depends on the identity (already pointed out in footnote 10 of Part II of this series) 


K[(ap + frp) -^2 (d _ y , (II Mb) ( II brj (f> 

P= 1 P= 0 V a GSd \i=l / \i=P+l 

holding for any d £ {1,2,3,...}, a pi b p £ R (p £ {1, 2,..., d}), where by convention we intend 

rui a <r(i) := n,= d+ i (j ) ;= i • 
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Summing up, Eq.s (B.9 


B.13) imply 




dr t s+ 2 * e 


d—p 


p ,o ( d -P)'P ! „ 6Si 


-I 1 B&,p( h) 

s-> _ r T~ 1 


(B-14) 


hez p ' 


T) ( Bcr.p(n) \ 

' s+^=2 V T / ■ 


EM 


To conclude and obtain Eq. (B.7), just note that all the integrals over r G (0,1) 


in Eq. (B.14) can be evaluated according to Eq. (3.33) to give the functions 


C Appendix. Absolute convergence and remain¬ 
der bounds for the series of subsection 3.7 


Let us recall that the expansion (3.79) for E u, ( > ' ) contains the series 


E 


neN d 


On the other hand, the right-hand side of Eq. (3.80) for E u, ^ < ' > contains 


(C.l) 


series m 


h G Z p , for p G {0, ...,d}. Let us consider one of these series; after removing the 
term with h = 0, which becomes singular at u = p + 1, we are left with 


E ^ir 1 

hezyh^o 


h) 


T 


(C.2) 


In the following subsections we prove the absolute convergence of the series (C.l) 


(C.2), for all u G C ; moreover, we derive remainder estimates for both series, which 


justify statements (3.85) (3.87). Throughout this appendix we adopt systematically 


particular, subsection A.l). 


the notations introduced in Appendix Al and the results obtained therein (see, in 


C.l The series (C.l); connections with the expansion (3.79) of E u ’(>) 


Let us define, for u G C and N G (0, +oo), 


77 


«>(>) ._ 


N 


E 

ngN d , |n|>iV 


^'“r( 


(C.3) 


If we can show that 77^ > ' ) < +cx) for some N , of course the series (C.l) converges 
absolutely. 
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Hereafter we derive quantitative estimates for 77^+ . Using the inequalities (A. 17), 
we infer from definition (C.3) that 


77 


«,(>) 

N 


< 


max(a SR “ 1 ,A* iu x ) 


7 r 


1 


E 

ngN d , |n|>iV 


n 


1—Sftu 


r 


t 


Kn —1 7T 2 T 


’ A 2 


n 


(C.4) 


now, using the result ( |A.13[ ) we obtain the following, for any a E (0,1): 

7 t 2 T Km —1 


v u,{» / niax(a R “ \A* U l ) {d) 

'^N + o d H N 


2 d yj-5Ru—1 


A 2 ’ 


, 1 —Km 


for either Km + 


l,N>2y/d or Km < 1, TV > 2v^ + -\/ —- : 

7 r V 2aT 


(C.5) 


Finally, let us point out that the remainder bound (3.85) for 77^/ > ' 1 follows easily 
from the above inequality noting that 


R 


«,(>) 

■N 


< 


K 


2|F(^)| Un 


u,(>) 


(C.6) 


C.2 The series (C.2); connections with the expansion (3.80) of E u, ^ < \ 
Let us define the following functions, for n E C, N E (0,+oo), p E {1,..., d } and 
a E S d : 


77 


N 


(RP) : = 

hgZP,|h|>JV 


*P u-p-i 


’5 CTiP (h) 
T 


(C.7) 


where T s and B a>p are as in Eq.s (3.33) and (3.81), respectively. 

Hereafter, we show that 77^ < ' ) (a, p) < +oo for some suitable N ; this implies the 
absolute convergence of the series (C.2). 

In the rest of this paragraph we derive explicit estimates for the functions 77^) ( ' < ' ) (a, p ), 
for all uE C, pE{ 1,..., d}, a E S d and some suitable TV. First of all, we note that, 
for all hGZ rf , from the definition of H CTjP ( h) in Eq. (3.81) it follows (compare with 
Eq. (fA24|)) 

a 2 |h| 2 ^|H CT;P (h)KA 2 |h| 2 ; (C.8) 

using these inequalities we easily infer (for all u E C) 


H ct , p ( h)V ^ max(a 5Rn “ 1 ,A 5R “- 1 ) |h 


I ffilL—1 


(C.9) 


Similarly to the derivation of Eq. (A.27), using the above bound along with Eq.s 
(A.3) (A.26) we obtain 


Ru-p-i 


B *A h ) 


T 


< 


max(a 3f?n p p x ) 

jRu-p-l 

7 2 


ijpUMt— p— 1 p 


/ p+1 —Km a 2 2 

V 2 ’ T |h| 


(C.10) 
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Inserting the previous results into the definition (C.7) of 77,[ t < ^(cr, n), we get 


max(a 5R “- p ' 1 ,A SR “-P- 1 ) ^ 


__ Ku—p—1 

T 2 ^ 


E 

heZP, |h|>JV 


( p+1 —Km a 2 2 




T 


(C.ll) 


using Eq. ( |A.5[ ) (with d replaced by n) , we conclude the following, for any a G (0,1): 

a 2 p+1 —Km 




max(fl s “ p 1 ,A 5R “ p x ) 


„ SRu-p-l 

1 2 


+?’ («, %; 


-,KM-p-l (C.12) 


1 /|Km —p—1|T 

for either Km + p+1, N > 2 y/p or Km > p+1, N > 2 y/p H— ' 


2a 


The above relation proves the finiteness of 77^^ (a,p) for all u G C, p G {l,...,d} 
and all a G Sd, which i mplie s the absolute convergence of the series (C.2); mo reover , 
the remainder bound (3.87) for R'^ 1 follows straightforwardly from Eq. (C.12), 
noting that 


R 


«,(<) 

N 


<^u T ^ * 

^ 


1 


E &a,p 

o 


77 


- (<) l 


a : p) . 


(C.13) 
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